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Abstract 

The process of multicomponent condensation is considered. The the- 
ory taking into account several channels of nucleation is constructed. The 
analytical approximate description of the whole condensation process is 
given. The specific phenomena of the secondary nucleation is described 
analytically. The possibility of the partial reverse condensation is dis- 
cussed. 
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1 Introduction 

Multicomponent nucleation is widely spread in nature and the task to describe 
it theoretically is rather actual. Creation of the classical theory of nucleation 
allows to investigate the case of multicomponent nucleation. The stationary rate 
of multicomponent nucleation was found in [2], solution of a stationary kinetic 
equation with appropriate boundary conditions was given in These results 
gave the opportunity to describe kinetics of the global evolution during the 
whole phase transition. The final period of evolution is known - it is the period 
of the recondensation (coalescence) and it is described by the Lifshic-Slczov 
theory (3]. But coalescence is mainly the redistribution of surplus substance 1 
and the real transition of practically all surplus substance occurs earlier. Before 
the coalescence, in fact, a phase transition takes place. To describe the phase 
transition it is necessary to determine the total number of droplets appeared in 
the process of nucleation. This quantity is the quasiintegral of evolution after 
the end of nucleation and before the coalescence. The characteristic times of 
coalescence ordinary strongly exceeds the characteristic times of experimental 

Certainly, during the coalescence a small part of the surplus substance is gradually going 
from the mother phase to the new phase. 
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observations and the characteristic time of a practically all surplus substance 
consumption exceeds the time of nucleation (i.e. the duration of the intensive 
appearance of droplets). Then the problem to determine the number of droplets 
is the central one and it will be solved here when many channels of nucleation 
a important. 

The problem to determine the global evolution until coalescence for multi- 
component case was investigated in [2], 0, [H], 0, JUj. Recent publications [8], 
HU repeat 2 the consideration made in [Jj], [7J, (TJJ- The mentioned publications 
gave the consideration of the most spread case of multicomponent condensation 
when the droplets of the sole composition will appear in the system. Then as 
it is noted in |5] , [7J , ^U] j the method of description of the nucleation can 
be reduced by renormalizations to the one component case. Therefore it was 
sufficient to give in 0, 0, E3) E! the corresponding references to the one 
component case ^B]- Nevertheless in 0, [H] in consideration of the binary con- 
densation the method of solution of one component situation was reproduced 
and unfortunately it was done with many errors. This requires to give here 
detailed analysis of this problem. 

Besides the consideration of nucleation with a sole composition of compo- 
nents in droplets the case of multiple compositions will be considered here. To 
see that the consideration based only on the sole composition isn't complete one 
we can imagine the following simple example. Let us suppose that the system 
contains two immiscible substances and both substances in a vapor phase are 
supersaturated. Then it is evident that the two separate processes of condensa- 
tion will take place 3 . As the result there will be droplets of the first substance 
and of the second substance 4 . So, it is evident that consideration of multi- 
component nucleation should contain the consideration of multiple channels of 
droplets formation. It will be done below. The effect of multiplicity of channels 
of nucleation is well known, it is deeply investigated by Ray [TT]. It is rather 
strange that in [5| (where ^I] is considered as the base for constructions) these 
conclusions of Ray are ignored. 

From the first pont of view it seems that this problem can be solved by 
the simple generalization of the case with the sole composition of substances 
in droplets. This illusion is wrong. The situation when there is absolutely no 
interaction between separate processes of nucleation of droplets with different 
compositions is exclusive. Normally there is such interaction and it has to be 
taken into account. This interaction forms the main content of multicomponent 
nucleation. 

To speak about kinetic effects of multicomponent nucleation one has to spec- 
ify external conditions. We shall consider the process of the metastable phase 
decay. In the initial moment of time we have a metastable mixture of gases. 

2 But contain some wrong statements, which will be corrected here. 

3 Certainly, these processes can interact by a heat release effects, but when these effects are 
small (it can be ensured by sufficient quantity of passive gas) one will observe two absolutely 
separate processes of the one component nucleation. 

4 Due to the entropy of mixing there can be some species of the other substance in droplets 
of the given dominating substance. 
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The stable phase is the liquid droplets phase, which will gradually appear in 
the system. Liquid phase droplets will consume metastable vapor molecules and 
the supersaturation of the mother phase is gradually falling. It leads to the ces- 
sation of appearance of new droplets and stops the process of nucleation. This 
is the qualitative picture of metastable phase decay, which will be described 
quantitatively. 

We consider the free molecular regime of vapor consumption. In |12| it is 
shown how renormalizations can reduce the situation with an arbitrary regime 
of vapor consumption to the case of the free molecular regime. Below some 
corresponding remarks will be given. 

The absence of thermal effects of condensation is supposed. Formally speak- 
ing this absence can be attained by sufficient quantity of a passive gas. The 
generalization on the thermal effects of nucleation is rather evident and can be 
done in a manner quite similar to [17| . 

We take the classical thermodynamic expression for the free energy of a 
critical embryo 5 . This requires that the critical embryo has to contain many 
molecules of every component. Only then we can use the approximation of a 
liquid solution in the embryo. Otherwise one can consider the molecules of a 
rare component as heterogeneous centers, but the requirement to have many 
molecules of all substances in the critical embryo remains 6 . Certainly, one can 
argue about the discrepancy between theoretical and experimental results in nu- 
cleation rates even in the one component homogeneous case, but it seems rather 
evident that this discrepancy can be eliminated by some fitting (normalizing) 
constants, which are very smooth functions of all ordinary thermodynamic pa- 
rameters of the systems including supersaturation and temperature. So, we 
can use the classical expression for the stationary nucleation rate as a brick 
for constructions made in this paper. It can be done because the variations in 
thermodynamic parameters essential for nucleation will be very small (details 
see later). 

The system is supposed to be homogeneous in space. The unit volume is 
considered, all energy like values are measured in units of a mean thermal energy. 

2 Channels of nucleation 

2.1 Asymptotic thermodynamic channels 

It is supposed that in the system there are I different substances, which can be 
found later in droplets of liquid. They will be called as components. Different 
components are marked by lower index. The molecular number densities are 
marked by rii, in the initial moment of time t they are known and equal to 
rii(O). The number of molecules of i-th component in the embryo is marked by 

5 Or some modified approach. As it can be seen below only some general features of the free 
energy (such as smoothness on the power of metastability of a mother phase) are required. 

6 The theory for heterogeneous decay is given in 1181 . 
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Vi, the total number of molecules - by 



* = 5>. 

i 

The embryo can be described by the set {vi} = V. Then the concentration & of 
component inside the embryo will be 7 

ti = vii$ ■ 

It can be done for every component and the vector £ = can be introduced 8 . 
Then one can speak about the molecular number density of a vapor of i-th 
component saturated over a liquid of concentration £ with a plane surface. It is 
marked by Uj e - According to standard definitions 

where n ioo = nj e (£i = 1) and /, are activity coefficients. Ordinary they satisfy 
conditions 9 

Ate -> i) - 1, -> i) -> o. 

In approximation of ideal solution = 1. 

The free energy of an embryo is marked by F. One can speak about F only 
for an embryo near the state of initial equilibrium. For the bulk liquid phase 
the free energy is given by Ft — — ^A/ij where the differences of chemical 
potentials are given by A/Xj = ln(nj/rij e ), which implies that the vapor phase 
can be regarded as an ideal gas 10 . Having extracted extensive variable i? one 
can rewrite F& as 

F b = -i?^^ln(n t /n !e ) 

z 

and extract the function 

B b = - 'y]€iln(ni/n ie ). 

i 

On the base of B}, one can formulate the necessary condition for the possibility 
of nucleation - function has to be negative at least at some £. 

The form Fb = i?B(,(£) leads to existence of the thermodynamic asymptotic 
channels of nucleation. These channels are straight lines in plane they are 
perpendicular to £-plane. The coordinates & c are extracted by conditions 

^ = o, 5 fl (fl t )-i;^(6-6c)fe-^c)>o . 

7 We don't take into account here the surface enrichment effects, ft will be done later. 
8 Here one concentration is absent due to y\ £j = 1. 

9 We shall use in this paper two sets of variables for the embryo description and 
{yi}. Sometimes ■& will be changed by the volume of the embryo. So, it is clear what variables 
will be fixed in partial derivatives. 

10 Certainly, virial corrections can be introduced here and this allows to describe many real 
situations. 
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These conditions can be satisfied at some £ c . Let the number of £ c be equal to 
lo- 

One can see that B b can be presented as 

B b = -"y]€iQn{ni/n ioo ) - ln& - In 

i 

The term — ln(^) ensures that 11 

dB b 

^-I^o = -oo, 
dB b 

This leads to the following properties of asymptotic behavior of F b : 

• (A) The asymptotic thermodynamic channel can not be at £i = 1 or = 
precisely. 

• (B) The previous property doesn't mean that B b goes to infinity, the 
equality xhix\ x ^,-\-o — > ensures that at £j = 0; 1 the function B b attains 
some finite values. Here it is supposed that all have regular behavior. 

The first property allows later to make some conclusions about the properties 
of nucleation channels. The second property establishes the possibility to include 
the unary nucleation into the multicomponent picture and to see the almost pure 
channels of nucleation, which will be important here. 

We shall denote coordinates of channels by £ as c . Consider now the binary 
case, we enumerate the channels to have 



as c 1 



< £j as c 2 < ••• < £,j 



Every channel has the zone of interaction. The boundaries of interaction 
zones are defined by relations 

d_B b= ^ 

The boundaries of these zones are straight lines in plane £, i9. They are perpen- 
dicular to £-plane. We have to get the coordinates of these lines. We shall denote 
these coordinates as £j as b j (the first index marks component, the last is the 
number, index "6" means "boundary"). The total number of these coordinates 
is lo — 1. They are distributed as 

£j os 6 1 ^ £j as b 2 ^ ••• ^ £i as b lo — l 

One can see that 

£j as c i—1 ^ Cj as b i—1 ^ £j as c it £,j as 6 i—1 ^ £,j as c i ^ £j as b i 



11 Under the regular behavior of /j. 
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The boundaries of the zone of interaction for the first channel are: the lower 
boundary is £j as & o = 0, the upper boundary is £j as b i- The boundaries of the 
zone of interaction for the last channel are: the upper boundary is £j as ; + i = 
1, the lower boundary is £j as / _i. The boundaries of the zone of interaction 
for other channels (number i) are: the upper boundary is £j a s b i, the lower 
boundary is ^ „« 6 i-i- 

One can determine the time of relaxation for every channel. It is necessary 
only to linearize the system of equations for the law of the embryos moving in 
plane vi near the coordinates of the channel. Then we have the ordinary dif- 
ferential equations with constant coefficients with known solution as relaxation 
exponents 12 . 

One has to note that for the two component case this picture is exact, but for 
multi component case the situation is more complex: in plane £ together with 
focuses we have saddle points. The analysis in this case here is the standard 
one for the theory of dynamic systems. 

The main ground for this analysis is the reduction of F into the product of 
intensive part B and extensive part i? ensures that the thermodynamic channels 
bottoms (when ■& is fixed) are straight lines. The same is true for the bounds 
of interaction zones - they are also straight lines in plane £, perpendicular to 
the £ plane. 

Also one has to note that even in the two component case there exist some 
illegal channels among channels of nucleation . They are situated in the region 
of the miscibility gap, which is defined by condition £?& > 0. 

One can give an interpretation of the function Bb as Afi - it is the mean 
value of chemical potential excesses 

3 

An interesting question, which will have important consequences is a ques- 
tion whether it is possible to see the asymptotic thermodynamic channel at 
some & w or at least at very small We shall these channels as the "pure" 
channels. The statements (A) and (B) show that it is quite possible but requires 
the specific behavior of activity coefficients One can extract the following 
properties of such channels. 

• (C) Dependence of fi over £ must be sharp. 

It means that for some £j the inequality \dfi/d£j\ S> 1 near the coordinate 
of the pure channel w has to be satisfied. 

• (D) The pure channels are narrow and deep 13 . 

The qualitative explanation is the following. The property \dfi/d£j\ 3> 1 
means that near the coordinate of channel bottom the absolute value of 

12 Here for the asymptotic channels the relaxation resembles the relaxation in the kinetic 
channel, which will be considered below. 

13 For the narrow channel the characteristic variations of £j are many times less than 1. 
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ft ordinary (not precise) attains big values. The existence of a channel to- 
gether with monotonous behavior of ln(^) at small concentrations means 
that fi has a minimum near the coordinate of a channel. The channel 
becomes deep. The non- linear behavior of ln(£j) at £j w ensures that 
the pure channel is narrow. 

• (E) Positions (coordinates) of channels are rather insensible to rij. 

Since fi has a deep minimum it means that slight variations of rii can not 
lead to essential variation of coordinates £, as c of a pure channel. 

The term "pure channel" doesn't mean that this channel corresponds to the 
absolutely pure condensation of one component. The entropy terms ln(£j) lead 
to existence of other components as species. This will be the main source of 
specific kinetic behavior discussed below. 

It seems now that the situation with the pure channel is rather rare in nature. 
Really the behavior of /» is required to be rather specific. But it is an illusion. 
Ordinary in the situation with immiscible components we have the mentioned 
behavior and the pure channels of nucleation. 



2.2 Real thermodynamic channels 

Now we shall include into consideration the surface term F s — Sa, where S is 
the surface of an embryo and a is the surface tension. The free energy of an 
embryo formation will be F = Ft, + F s . This form is the so-called capillary 
approximation. All other correction terms are neglected here. The term F s is 
asymptotically negligible. 

We shall try now to repeat all constructions for F instead of Ft,. Really it is 
possible to do as far as F can be presented in the following form 

f = +nB(s) + n 2/3 , 

where fi = (Sa) 3 / 2 , then 



B(0 = 



^3/2(^6^1/2^ 



Here v is the mean value of the volume per one molecule in a liquid phase defined 
as 

j 

and vi j is the volume occupied by a molecule of j-th component in a liquid 
phase. Then 

F c = fi 2 / 3 /3 

(index c marks the critical embryo) and the ordinary Gibbs rule takes place. 

Later we suppose that there arc no irregularities in behavior of % and 
a as functions of £. 
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Now we can see the maximum over the extensive variable Q and again we see 
the set of channels defined by the extremums of function B. We can reproduce 
here all arguments and constructions as we have done for F^. Again we observe 
properties (A)-(E). Certainly, the numerical values of coordinates of channels 
will be another. We shall call these channels as "thermodynamic channels of 
nucleation" . 

The difference between consideration of the real thermodynamic channel and 
the asymptotic thermodynamic channel is only in the process of relaxation. For 
the asymptotic channel one can not observe the saddle point in thermodynamic 
channel. So, there is no use to speak about the quasi equilibrium along i? in the 
precritical region. For the real thermodynamic channel one can speak about the 
quasi equilibrium. Then one can estimate the time of relaxation t r over £ in a 
way it had been done by Zeldovitch for the near critical region and come to 

Also one has to add to this time the time of regular motion to the viscinity of 
the "bottom" of the channel, i.e. from initial values to the region |£j — £j c | ~ 
v i y 2 \d 2 B/dti t f\~ 1 / 2 2 1 / 2 . Since B is arbitrary, we prefer not to give the final 
expressions here, because it can be done only on the base of the behavior of B 
after the evident integration. 

Here it is necessary to outline the meaning of v. One has to keep in mind 
the Gibbs absorption equation or the Gibbs-Duhem relations. In application to 
the embryo description it means the phenomenological rule, which states that 
having come to v and £ one can simply neglect all derivatives of v over £j . The 
same has to be done with derivatives of a over & , which is known as the so-called 
Wilhcmsky-Renninger effect |13|. This rule, which is very useful in practice, can 
be shown on the base of Gibbs absorption equation. 

In transition from F& to F the surface tension have appeared. Surface tension 
can have rather arbitrary behavior on concentration. Can it violate the previous 
conclusions? The typical behavior of a as function of £ is the following: as 
long as possible there exists an enrichment of a surface layer and the surface 
tension doesn't variate essentially. Certainly, the entropy of mixture has to 
be contained in a. Only when there is practically no molecules of component 
with a small a, then the value of a begins to grow and attains for the pure 
substance the value corresponding to a for this pure substance. It means that 
the pure thermodynamic channel of nucleation for component with maximal a 
is practically impossible 14 . 

This picture of the surface tension behavior has one important qualitative 
consequence. Precisely speaking the asymptotic thermodynamic channels have 
nothing in common with precise thermodynamic channels and are useless. Cer- 
tainly, for the ideal solution (all fi = 1) the additional entropy of mixture is 

(jlsurf I (jlsurf ~t~ sur/)) ln(fi sur / / [71 sur / -\- Tli surf))~\~ 

surf I (jlsurf ~1~ sur/)) ln(f^ surf/ij^surf ~t~ T^i sur/)) 
14 Certainly, the sharp dependence of a can be compensated by peculiar behavior of /j. 
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(nbuik/ (nbuik + rii bulk)) ln(nfe„ife/ (n bu ik + rn bulk))— 
(«» bulk/ {nbuik + n» bulk)) ln(rii («6«(fc + «i buifc)) 

where sur / is the density of given component in the surface layer, rii bulk - 
the same value for the bulk phase, n sur f - the density of all other components 
in a surface layer, nbuik - the same value for the bulk phase 15 . Then the sum 
over all components has to be taken. We see that the additional entropy of 
mixture at moderate £j cannot be a sharp function of concentration. Under the 
rough generalization one can speak about the flat region - the region of the flat 
behavior of a as function of and about the sharp region - the region of the 
sharp dependence of a on £j. In the sharp region the thermodynamic channel of 
nucleation can be hardly situated, then this region isn't important. In the flat 
region the derivative of surface tension is small and it means that the position 
of the thermodynamic channel can not be essentially changed by inclusion of 
the surface term. So, one can use the asymptotic (bulk) thermodynamic chan- 
nels as initial approximation to find the coordinates of the real thermodynamic 
channels. 

One has to show inclusion of the picture obtained in consideration of Fb into 
the picture in consideration of F when $ goes to infinity. Precisely speaking 
the roots of B don't coincide with the roots of B b . This is because the roots 
of B are coordinates of channels when the extensive variable fl 3 ^ 2 is fixed and 
the roots of Bb are coordinates of channels when $ is fixed. Certainly, these 
coordinates will be different even when the channels would be the same ones 16 . 
To compare channels one has to go from to i? according to 

Then the necessary inclusion will be evident. 
2.3 Kinetic channels of nucleation 

Now we shall turn to the definition of kinetic channels of nucleation. Really, the 
mentioned thermodynamic channels of nucleation arc important for transition 
over the activation barrier, but when the droplets are moving along the slope of 
the free energy profile then kinetic terms are extremely important. At the sizes 
strongly greater than the critical sizes one can not speak about the absolute 
extremums and saddle points, which are independent from kinetic factors. Now 
we shall investigate embryos with Vi c more than three-four times greater than 
the corresponding critical values. 

Now we have to investigate kinetic equations and to determine their station- 
ary points. The simplification here is the possibility to use only the regular laws 

15 Here we suppose that the function v(£) is very smooth. 

16 When H — * oo or ■& — > oo then \F — Fi,\ -C F, but this doesn't mean the similarity of 
coordinates of channels. 
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of growth 17 for embryos (since they are supercritical embryos here they will be 
called as droplets). 

The regular laws of motion 

— = W + - Wr 

dt 1 1 ' 

where and W[~ are direct and inverse flows on the embryo (the coefficients 
of condensation are included in these flows) lead to expression for 

^^-f-tf- 1 V^+ ^ r 1 — 
c^^ M ^ dt dt 

3 

So, the coordinate of a kinetic channel £j k is determined 18 by 

dvj dt/i _ 

« dt + dt ~ u ■ 

j 

It is necessary to stress that the last equations will be equations on £, the 
dependence on z9 is absent here. Really, one can write for the following 
formulas from a simple gas kinetic theory (it is supposed that gas is an ideal 
one and the regime of substance exchange is the free molecular one) 

= ^Vi tSrucn 

Here Vi t is the mean thermal velocity of a molecule of component i, on is 
condensation coefficient 19 . The surface square of an embryo S is the same as 
in the expression for the free energy, it differs only in microscopic corrections. 
The kinetic coefficient W~ can be obtained on the base of by the detailed 
balance relation 

Wr(0, Vi) = W+(v, Vi - 1) exp(-F(j7, Vi - 1)) exp(F(i?, Vi )) 

In the supercritical region near the bottom of the channel the function exp(— F(y, Vi — 
1)) exp(F(i7, Vi)) doesn't depend on ■& when V = J2 i ViV t is going to infinity and 
one can easily see that d^i/dt can be presented as i? multiplied by some given 
function of £. It allows to speak about the line as a limit of kinetic channel in a 
plane #, £. Now we shall determine this line concretely. Instead of the previous 
relation one can write 

Wr(v, v.) = W+(v, Vi - 1) exv(dF(y)/dvi) 

17 This requires a special justification for the supercritical embryos - the embryos, which 
contain much more molecules of every component than the critical ones. This justification 
can be given analytically. 

18 Here and later the lower index k denotes the stationary characteristics in the kinetic 
channel. 

19 It can be included into the mean thermal velocities Vi f 
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Calculation of dF(y)jdvi leads to 



dFW/Ovi = -ln( n J— ,) + (for 1 ' 2 ) 2 ' 3 2 (^v^-^avj + G.D. 

Here G.D. marks the terms, which according to the Gibbs-Duhem relation have 
to be cancelled 20 . 

The Gibbs-Duhem equation for this case can be written as 

3' 3 J ' j' J 3 

for arbitrary i. This Gibbs-Duhem relation isn't the standard one because it 
contains the term, which has a low dimension. This equation appears as a 
combination of the Gibbs-Duhem equation for the bulk phase and the Gibbs- 
Duhem equation for the surface layer. One has to keep in mind that due to the 
surface enrichment the concentration & in a bulk solution of an embryo can not 
be precisely defined as Vi/d. This note leads to the last relation. 

One can see that the term (67T 1 / 2 ) 2 / 3 iQ^ Vtv^^^aVj is negligible in com- 
parison with \n(rij/(nj oo€jfj(€))) when the volume of an embryo ViVi is big 
enough. So, one can neglect the correction term and come to 

r> ■ 

dF<p)/foi = -H f—r) 

This relation is typical for asymptotic kinetic channels of nucleation. Here one 
can already speak about the straight channels of nucleation because the rhs 
doesn't contain extensive terms and dvj/dt can be presented as 5 multiplied by 
a function of intensive terms. 

The Gibbs-Duhem relation becomes here also more standard and can be 
written as 21 



EE< 







This form is absolutely typical for the bulk solution with standard definition of 
concentration £j as Vi/-d. 

For asymptotic kinetic channel the relaxation to the stationary value of con- 
centration is practically evident. The absence of curvature here is very impor- 
tant in contrast to the thermodynamic channels in the under-critical or near- 
critical regions. The matter is that the variable •& in the super-critical region has 

20 One can argue whether the Gibbs absorption relation can be applied to a growing embryo, 
but we have to adopt this equation in order to have concrete calculable results. 

21 Here appears a special question whether it is possible to consider Vi as and to take 
i?c9[£j ^{n-j /njac^j fj)]/d^j instead of d[uj ln(n j /n joo£j fj)]/d£j m two previous Gibbs-Duhem 
relations. Since the channels of nucleation are rather narrow (£j — £j c ~ i? -1 / 2 ) this can lead 
only to correction terms with negligible order in powers of ■&. 
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to be regarded as the fast variable. This lies in contradiction with consideration 
of the near-critical region (where are the fast variables and t? is the slow one). 

Consider the near-critical region. Having calculated dF/dVt and dB/d^i 
one can easily see that if there are no small parameters in , fi and their 
derivatives, then the relaxation to the bottom of thermodynamic channel takes 
place and the characteristic time of this relaxation is estimated as ~ -d 1 / 2 jW^ . 
The time to overcome the barrier in the near-critical region has the estimate 
~ $ 2 / 3 /W^~ and one can see that there exist the hierarchy between these char- 
acteristic times, which provides that along the stable variable ~ — & c ) there 
is the quasi equilibrium. 

The opposite situation takes place in the super critical region. As we shall 
see later the variable i9 grows very rapidly, this property will be called as the 
avalanche consumption of metastable phase (see \'2'2\ ). When the channel is 
straight one then there is no difference whether i9 is the fast variable or not and 
one can see such a relaxation. But when the channel is curved there can be 
no relaxation to the bottom (the term "bottom" at the slope becomes rather 
relative). So, here the asymptotic disappearance of curvature plays a very im- 
portant role. Really, the velocity 22 dfi/dt allows the estimate 23 

dti/dt ~ tf 2/3 -> oo, 

- & k)/dt ~ ti^dvi/dt ~ ?T 1/3 -> 0. 

Fortunately ■& ~ t 3 and 

d& - 6 fc ) 



dt 



dt ~ / dt/t ~ Int 



So, relaxation really takes place, but it isn't the exponential, but the power 
relaxation 24 . Here the straightness of the channel plays the principal role. 

We shall enumerate kinetic channels by the number of thermodynamic chan- 
nels, which this kinetic channel will be attained from (So, the kinetic channel 
can have several numbers or no number at all. Kinetic channel can be described 
by one of its possible numbers.). 

One can not investigate in details the transition from the critical region to 
the region of sizes, which strongly exceed the critical sizes. One of effective 
approximations is the 5-function approximation suggested in 0. Let some sta- 
tionary value of concentration £j (it is one dimensional variable for the binary 
system) be & + . The characteristic width A£$ along £j allows the estimate 

A& ~ zT 1 / 2 . 

The last value is relatively small. So, in the law of growth one can put & w £j + . 
Then 

^ „ ,92/3 

dt 

22 When the embryo is moving near the stationary value, i.e. the concentration alterations 
are small. 

23 Here and later the lower index k denotes the equilibrium value in kinetic channel. 
24 All details can be seen in 6 , Appendix 1 and Appendix 2. 
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which leads to $ ~ t 3 . The characteristic deviation in plane vi can be estimated 
as (£» ~ & +)& an d has a characteristic value i? 1 / 2 . Having rewritten the last 
estimate as a function of time we get (£, — & <~ t 3/<2 . We see that this 
width grows faster than the width of the diffusion blurring ~ i 1 / 2 . Then one 
can consider the pure diffusion over £. The solution of this problem is known. 
The distribution n over #, £ with 5-like initial conditions 25 will be presented as 

n ~ n# n$ 4 , 

i 

where n# ~ J d^n is simply a translation of initial conditions (i.e. the boundary 
conditions at small 26 d) 27 and ng 4 approximately satisfies the diffusion equation 

dntJdt~ Di&nzJdg 

with characteristic diffusion coefficient 28 Di and looks like the Green function 
of diffusion on the flat energy profile 

G ~ exp(-(& - & +) 2 /(4D i t))/(4*D i t) 1 ' 2 

Approximation Dj w const as the function of size, which was used to get the 
previous relation isn't too suitable. Now we shall extract coordinates with a 
real physical meaning. 

Having noticed that i? <~ t 3 , <9/9i <~ 9/di? 1 / 3 we can write instead the previ- 
ous diffusion equation the following equation 

dntJd#~ D 0i d 2 nii /d^ . 

Here D i is the new diffusion coefficient for evolution in £ plane. The really 
significant variables are namely these variables. A solution of the last equation 
has a good approximation also in the form of the Green function 

G ~ cxp(-(& - & + ) 2 /(4A>^))/(4ttA>^) 1/2 • 

The diffusion over the plane surface can not lead to the approaching of £j 
to £j fc. To see this approach one has to take into consideration the terms 
corresponding to the regular growth. 

For £j near £j k one can restrict the regular law of growth only by linear terms 
in (£; — £j fc) series. For diffusion terms one can take the ordinary coefficient 
of diffusion. So, we observe the evolution in a square potential. Then one can 
prove that it is possible to use the Fokkcr-Planck approximation. 

The Fokker-Planck equation in multi dimensional harmonic potential has a 
well known solution. Multidimensional square potential in £ is invariant over 

25 Because the near-critical region has the small relative sizes. 
26 But still greater than critical. 

27 In the stationary conditions it is the stationary rate of nucleation divided by the regular 
law of growth dd/dt 

28 All W + arc supposed here to have one and the same order. 
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Lorenz transformation, which leads to the separation of variables. We shall 
denote the new variables also by £j here. 

Since Lorenz transformation is a linear one, it can not lead to any terms 
except linear ones along Now we have (I — 1) separate diffusion problems 
with known solutions (the Green function is known, the eigenfunctions and 
eigenvalues of harmonic oscilator are known, one can act in both approaches) . 

2.4 The coordinates of channels 

Now we shall speak mainly about asymptotic kinetic channels and forget about 
the curved kinetic channels. The reasons will be the following: 

• When the channel is curved it means that the current characteristics aren't 
too far from the critical ones. Then the the derivatives dF/dvi are small 
ones and one has to take into account the diffusion. The interesting situ- 
ation is far from the critical region where the channel is straight. 

• The characteristic size of relaxation to the stationary value of concen- 
tration strongly exceeds the distance where the curvature of a channel is 



• Later namely the characteristic sizes corresponding to the asymptotic 
channels will be important for construction of the global kinetics of the 



We shall mark the coordinates of asymptotic kinetic channels also by the 
subscript k. We shall call the asymptotic kinetic channels simply as the kinetic 
channels. 

To determine this coordinate of equilibrium in the kinetic channel one has 
to write the laws of regular growth 



at the supercritical asymptotics. 

The surface square S is given by S — (67r 1 / 2 ) 2 / 3 u 2 / 3 (^)i? 2 / 3 , which leads to 



essential. 



process. 



dv l 
dt 



d£j _ (6^/2)2/3 



v 2/3 [aiVi t (rii - Hi 



dt 40V3 



'l oo 







Here one can extract the function 




Equations 



Ki = 



for all % together with 
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form the system of nonlinear equations for coordinates 29 £j k- This system has at 
least one solution when the vapor mixture is metastable. It can be derived from 
the boundary behavior. But sometimes this system can have several solutions 
(this corresponds to several channels). 

Near one can rewrite the previous equation as 

(6^ 1/2 ) 2/ y 1/3 ^(a , ,-2 /3| 



Equation for $ will be 

M _ (67T 1 / 2 ) 2 / 3 
~dt ~ 4 



or 

where 

( fi l/2\2/3 

One can suggest the simple approximation when & is near & & for all i (here 
some given channel is already chosen) 

This equation can be easily integrated 

,9 1 / 3 (t) = (t-to)P(er fe )/3 

where to is the time of appearance of a droplet in a supercritical region 30 . 
It allows to rewrite the equation for d^i / dt as 

d£i _ (eTr 1 / 2 ) 2 / 3 1 ^.(^ ^ ^ 



* 4 (t-to)P(^ k )/S d£i 

with solution 



(& - & fe) = const(t - t ) 

and unknown const where 

(6.v 2 ) 2 / 3 i dm) 2/3| 



29 It is necessary that the square form ^\ - (d 2 F / d(,id^j)\k{(,i — €,i k){(,j ~ £j fe) wm be 
positively defined. Later only these £j ^ will be considered. 
30 Initial size can be put to zero. 
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If one can spread the last equation until S,i = £,i c (index c marks the critical 
value) one can very approximately get 



const = (& c - & fe )/(i c - t ) L 



where & c is the concentration of the critical embryo (one has to decide, which 
critical point has to be chosen), t c is the time when according to the regular law 
of growth the embryo attains the critical value of #. 

One has to stress that there appear at least three problems, which have to 
be solved 

• 1. how one can get the critical values of •& and £, i.e. d c , & c ; 

• 2. how one can get the asymptotic values of and £; 

• 3. how one can decide, which asymptotic channel corresponds to the given 
critical point (i.e. to the given thermodynamic channel). 

The first two problems are purely algebraic ones. One can only note that to 
solve the first problem one has to use the Kelvin relation instead of direct search 
of minimum of function B because of deviation of the bulk concentration from 
the integral definition, which occurs due to the Gibbs absorption isotherm. 

As an initial approximation for solution of the first problem one can take 
the roots of the function dB/d^i, where it is necessary to forbid the differenti- 
ation of a. Then one can see the regions of interaction for all channels having 
calculated d 2 B / 'd&d^j and decided whether the suspicious root is the bottom of 
the channel, the top of the hill or the saddle point. Then dB/d^i = leads to 



which is the known Kelvin relation. 

Function Bb can be found from B when Vi have approximately the same 
values and a is independent on the concentration, sometimes the search of the 
roots dBb/d^i = is more simple than the search of the roots dBb/d^i = and 
can be regarded as some initial approximation. When all vi have approximately 
the same values then asymptotic thermodynamic channels lies near the real 
thermodynamic channels. 

To fix the positions of asymptotic kinetic channels of nucleation sometimes 
one can use the approximation of strong supersaturations. This approximation 
means that in Ki the value rii oo£i/i(£) is neglected in comparison with m. Then 



v 



dAu 




This leads to 



Vi Vi V 




3 
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and together with £i = 1, (which now is trivial) we have the system of linear 
equations with known solution 

_ OLiVj tUj 

si v-^ 

L, a j V j tflj 

Here the solution is unique. 

Another promising approximation is approximation of ideal solution. In this 
approximation one can formally split the system into two components: the first 
component and all other components. Then one can get kinetic coefficients for 
these two components and solve the two-component case where the equation 
of the concentration of the first component will be the following (obtained by 
D.Stauffer in 0) 

aivi t (l - - ni 00C1) = "2^2 t£,i(n 2 - n 2 oo(l - £i)) 

(it has two solutions: one for the bottom (stable stationary value) of the chan- 
nel and one for the top (unstable stationary value), also it is necessary to check 
boundary values at £i = 0; 1 as possible channels). The problem is that n 2 oo 
is unknown because the composition of the second complex component is un- 
known, but now for this component the number of substances is reduced and we 
can continue this procedure having chosen in the second component two com- 
posing components, etc. As a result we shall get the equation for concentration. 
This completes the method of approximation of ideal mixture. 
In the general case one can speak about invariant 

aiVi t (rii - m ooiifi)/U = inv 

and write equation on this value. This way was described in |18| . 

The third problem is rather hard to solve but one can note some essential 
properties: 

• Ordinary one can use the regular velocities of droplets growth and ig- 
nore fluctuation terms. Only under a very specific behavior of properties 
the fluctuation terms are essential. This can take place only when there 
explicitly exist some small parameters. 

• Rather often it is possible to act in a following simple manner: It is suffi- 
cient to get the regions of interactions for kinetic channels and then to see 
to the region of what kinetic channel belongs the given thermodynamic 
channel. This way is very effective for practical needs. The probability to 
make an error is very low. 

• In any case it is impossible that droplets can cross the miscibility gap. 
Results presented in H are wrong. Really, in the miscibility gap the slope 
is positive (with increase of ■& the value of free energy also increases) . So 
here the kinetic channels will be thermodynamic ones (if we don't take into 
account the disappearance of the embryo), which have no saddle points. 
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3 Standard nucleation 



Nucleation is a period of a rather intensive formation of droplets. To begin the 
investigation of the nucleation process one can note that the number of ther- 
modynamic channels have nothing in common with the number of compositions 
of solution in droplets, i.e. with the number of kinetic channels. The illusion 
of such connection was presented in [§], but there are no reasons to justify this 
illusion. 

According to the ordinary adopted classification one can extract in the con- 
densation the first period - the period of intensive formation of droplets and the 
second period - the period after the nucleation but before the coalescence. But 
it is necessary to stress that even at the first period the variations of the vapor 
densities n, or of activities Q are rather important. This is in contradiction with 
the statement made in [9] (page 391) where it is noted that rij during the period 
of nucleation are constant. Under such assumption it is impossible to get any 
spectrum of the droplets sizes 31 different from the constant amplitude. Nev- 
ertheless in 9 some attempts to get the spectrum of sizes were made without 
any citations of the already presented theory in p], JU|, [7]. They follow 
instructions of [5], E3j EH > El to consider the one dimensional equation instead 
of the system of condensation equations, which also were presented in [5], [TU]. 

As it was noticed earlier & £j It means that the quantity of the 
molecules of component i in droplets (i.e. <jfj) satisfies equality 

9i/£i k = 9j/€j k 

for two different components. This takes place when the characteristic size of 
droplets during the nucleation period (the first period) strongly exceeds the crit- 
ical size and one can assume that the evolution is mainly determined by rather 
big droplets, which have £j ~ £,i k- The process of relaxation was described 
above. 

Namely this relation is the base for reduction of multi-component case to the 
one component condensation. The relaxation processes haven't been mentioned 
in [5] at all. Really, having marked the initial densities of vapor as rii(t = 0) 
and the current densities as rii(t) (t is time) one can easily see that 

(n;(0) - n<(t))/6 k = K(0) - n,-(f))/& k 

for two different components i and j. This forms the base for reduction to the 
one component case. 

We see that this reduction isn't the consequence of some functional approx- 
imations used in [§], ^U] and repeated in [§]. So, it is more than symbolic 
that to get the reduction to the one dimensional case [5] also uses the func- 
tional exponential approximation and approximation of a regular growth (they 

31 Here and later the size of the droplet will be the characteristic, which grows with the 
velocity independent of this characteristic. For the free molecular regime such characteristic 
1 /3 

will be f ' . The concentrations here arc fixed. 
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are also from [5], (TJU, which isn't necessary here and is the little defect of 
narration 32 of [5]. |rj]. [TU]. 

One can directly show that the droplets of big sizes are the main consumers 
of vapor phase. One has specially stress attention on this property despite it 
was considered in JU] because in [5] in was announced (p. 391) that "this hardly 
influences the results of our theory but significantly simplifies the mathematical 
expressions". The cited conclusion is wrong 33 . This property can be shown 
by the set of simple qualitative estimates. We omit here these estimates be- 
cause this property will be directly seen from the final results. Also it can be 
directly shown that the quasi-stationary approximation for the droplets, which 
sizes i? two or three times exceed the critical size # c really takes place (the con- 
centrations here are arbitrary, but in reality the quasi-stationary distribution is 
essential only for concentrations near the stationary value in the given channel) . 

Both these properties are based on the possibility of the thermodynamic 
description of the critical embryo. In 9 the possibility of analytical derivation 
of these properties hasn't been mentioned at all. 

On the base of the last statements one can write the balance equation in the 
following form 

7*(0) =n,-(t)+ft(t) 

where gi (t) is the number of molecules in the liquid phase and it can be presented 
as t 

g t (t)= [ dt'piityifatyt 
Jo 

Here Vi(t,t') is the number of molecules inside the droplets born at t', pi(t') 
is the intensity of formation of droplets. Our next task is to determine this 
intensity. 

The intensity of the droplets formation depends on the height of the acti- 
vation barrier of nucleation. Since there are several saddle points it seems that 
there are several possibilities for embryos to become the supercritical objects of 
the liquid phase. Ordinary, the leading channel is the unique one. The intensity 
of formation is proportional to exp(— F c j) (here F c j is the height of activation 
barrier in the j-th thermodynamic channel). It follows from the standard clas- 
sical consideration (see POD in stationary conditions because the nucleation rate 
can be calculated in stationary approximation. Here it is sufficient to keep only 
the leading term and not to think about the preexponential factor, which is the 
matter of numerous discussions 0, 0- It is important that the normalizing 
factor of the equilibrium distribution for every thermodynamic channel is one 
and the same since all channels go from the origin of coordinates. 

One can choose the channel with the lowest F c j. Let it be the first one. 
If there is another channel with | F c j — F c i \ < 1 the theory presented below 

32 It this derivation it is supposed that the equilibrium concentrations do not change too 
fast, which has to be proven (and can be done). 

33 It isn't clear how after this conclusion in |5] it was possible to consider only the droplets 
of big sizes without the logical destruction of the theory. 
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has to be slightly modified 34 . But the probability of such coincidence is very 
low since F c j ^> 1. All these estimates have to be done at initial values rii — 
Ui{t = 0) when there is no depletion. At first we shall investigate the period of 
intensive formation of droplets in the first channel. This period doesn't cover 
all nucleation possibilities but as we shall see later namely this period will be 
the main in nucleation. 

From the first point of view it seems that the variations of n, can lead to the 
variations of AFi c and sometimes the channel, which was the first one can be 
changed by some other channel, which becomes the most energetically profitable 
one. Fortunately the leading channel remains the main one at least until the 
end of intensive nucleation in this channel 35 . 

The droplets consumes vapor and this effect strongly diminishes the intensity 
of nucleation rate. To see this one can simply calculate the derivatives of the 
free energy of the critical embryo 36 

dF c 1 

"3 — = v i c 

dn t 

The last relation can be seen directly from the explicit expression for the free 
energy F(v) if we notice that in 

dF c _ dF ^ dF dvj 
drii drii ^-r 1 dvj drii 

3 J 

the term V M^-S^- vanishes as far as M^- — at the point of extremum (i.e. at 
the saddle point). This equality is a concrete realization of nucleation theorem 
in the multi-component nucleation (see [Hj, for the first derivation in the 
application to the multicomponent condensation kinetics, but generally this 
result was mentioned in [T%|). 

This relation plays the leading role in further analysis. To apply the ther- 
modynamic description of the critical embryo it is necessary to have ?? ^> 1. 
Without the thermodynamic description it is nothing to do in the quantitative 
theory of nucleation. To use the approximation of solution in the bulk of the 
embryo it is necessary to be c 1 . Otherwise we have to construct the model 
of liquid solution in the force field of several molecules of rare components. This 
problem is very hard to solve. Then it is treated as heterogeneous condensation 
and the molecules of rare components are considered as heterogeneous centers. 
The theory of heterogeneous condensation has been given in |15| . So, here we 
have to suppose that all i>i c 3> 1 (but this doesn't mean that & essentially 
differs from zero for the supercritical embryos). 

On the base of —riidF c / 'drii = i/j 3> 1 one can easily see that even the 
small fall of the molecular number density leads to the interruption of intensive 

34 We have simply to add this intensity to the intensity in the first channel. Functional 
dependencies on the vapor densities will be approximately the same. One has to take into 
account the intensity with the smallest Vi c (see later). 

35 The theory for this exclusion will be presented later. 

3e The thermal effects are supposed to be absent as it has been noted earlier. 
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(in comparison with initial) formation of droplets. It is the base for further 
approximations . 

One can see that this point of view radically differs from opinion presented 
in |5] where inequalities like dF c /drii ^> 1 are considered (see [5] eq. (23) ) as 
some special artificial external restrictions. 

Since riidFj c /drii 3> 1 then the main dependence of the stationary intensity 
of formation J s , of the droplets in the i-th thermodynamic channel is given by 

J s !(n(0)) = J s ttnfi)) exp(-Fi c (n(0)) + F 1 c (n(t))) (1) 
This formula is for the first channel. It is valid when 

K(*) - ni(0)| < ni(0)/ui c i| n4=rii (o) (2) 

This region of n.; corresponds to the intensive formation of the droplets through 
the first channel. The last estimate restricts variation of n.j and the process 
of formation takes place in time. One can imagine that a long tail with a 
small intensity can be more important than the short initial stage with a high 
intensity. Then the last estimate fails. Fortunately it is easy to see that every 
droplet grows with a growing intensity in time. Then d 2 rn/dt 2 < and the 
estimate (jSJ really corresponds to the region essential for nucleation. 
One can see one very important feature for rii satisfying J3J): 

• For m satisfying J2J all values of £j c i are approximately constant 37 (the 
relative variations of £j c i and 1 — & c l are small). For all other thermo- 
dynamic channels the relative variations of equilibrium concentrations are 
small if the channel is the "clear" one 38 . 

The analogous property takes place also for kinetic channels 

• For rii satisfying @ all values of & i are approximately constant (the 
relative variations of k l an d 1 — £j k l are small). For all other kinetic 
channels the relative variations of equilibrium concentrations are small if 
the channel is the "clear" one 39 . 

These properties can be justified analytically. 

The estimate (01 also allows to make decomposition in the argument of 
exp(— Fi c (n;(0)) + F\ c (ni(t))) in Tailor series. It is necessary to decide what 
number of derivatives has to be taken to ensure the high accuracy. In [5] the 
parameters of analogous decompositions were chosen in a wrong way. Really, in 
P| the correction term has an order proportional to the first derivative dF c /drii, 
the second derivatives were ignored. It is wrong and has a hidden supposition 

37 The last index is the number of a channel. 

38 The "clear channel" is the channel with relatively high sides (more than several thermal 
units). 

39 Here the sides of the channel are determined i9 being fixed. 
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that the second derivative have the same order as the first derivative. This 
supposition is too strong. Instead of this one has simply to require that 

dF 

-Fi c (n,(0)) + F 1 c (m(t)) « -^1^=^(0) " «i(0)) 

The last approximation has to be directly checked for every type of multicom- 
ponent mixture. It can be shown on the base of the evident inequalities 

-> -> dF 

I ~ Fi c(«i(0)) + Fi c(ni(t))\ > I^^I^^^W-a^O))! 

i % 

and 

\-F x c (n i (0))+F 1 c (ni(t))| < | £ -^U^oxi-i/"^ k=Bi(0)) M*)-»<( ))l 

There is no special parameter of decomposition, the considerations made 
in O-symbols presented in [H] are irrelevant because in reality the values of 
dF\ c /dm are big, but still finite. Nevertheless it is possible to give some con- 
vincing arguments to see that the last approximation is valid. We shall show 
that the terms with second derivatives are small. For one component condensa- 
tion the required property is obvious and can be checked by explicit formulas. 
In the multicomponent nucleation since there are generally unknown functions 
fi we can not get the general results. To estimate the second derivatives we 
have to consider the values (dvi c /dnj)/vi c . They take moderate values, at 
least they are smaller 40 than 2/3 c . To see this property one can use variables 
i9, £ and if all fi are directly independent on rij one can see that [dvi c /drij)/ui c 
has no small or big parameters. So, (d^ c /drij)/vi c takes moderate values less 
than the values c . 

According to J2J) in the law of growth one can take n, = n,i(t = 0) and get 



dv) IZ 1 



dt 12 

The values of concentrations £ have to be taken as the equilibrium values in 
kinetic channel corresponding to the first thermodynamic channel. Here we 
shall take more strong definition of the supercritical embryos. We shall call the 
supercritical embryos as those embryos, which are already passed the critical 
point and have already approached to the kinetic channel. 

One has to check that even with a new definition the statement that the 
supercritical droplets are the main consumers of vapor remain valid. Under the 
conditions of applicability of thermodynamic description of the critical embryo 
this statement can be proven analytically with some reasonable assumptions. 
Also the second statement about the quasistationary value of the nucleation 
rate during the nucleation period has to be proved for a new definition of the 

40 Here it is supposed that all rij have one and the same order. 
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supercritical droplets. It also can be done analytically in conditions of thermo- 
dynamic description of the critical embryo. 

One can state that the same property as for thermodynamic channels can 
be seen for kinetic channels. Namely for variations of n, satisfying J2J the 
equilibrium concentrations of kinetic channel can be considered as constant 
ones. This property is valid for the "clear" kinetic channels when that the 
derivatives d^i/dt have no small parameters in the vicinity of the equilibrium 
values in kinetic channel. 

For densities rij satisfying to J2J one can take also £ as constant values £1 & 
and get the rhs as a constant value. Then it is possible to integrate the last 
equation and get 

vl /3 = Mt-t') 

where 

Ai = a i ^(67r 1 /2)2/3^2/3^-2/3 u . t(n . (Q) _ n . ^ i fe/ .(^ fe) ) 

Here it is taken into account the initial condition i>(t = t') = 
It is easy to see that 

Af 

£,i 1 k £,j 1 k 

for different components. 

The equation for gi will be the following 

t 

dt'J{t')Af(t-t') 3 

Here J(t') is the rate of nucleation. 

For J one can take according to the second statement the following expres- 
sion 

J ~ Js 1 

where J s ± is the stationary rate of nucleation through the first thermodynamic 
channel. 

Then approximation Q leads to 

J = J*exp(-^—j-^\ ni = ni ( )(ni(t) - n;(0))) 

i % 

where 

J* = J s i(n = n(t = 0)) 

Then 

/•t jp 
9J= I *'J»exp(-^^|„ i= „ i(0) (n l (t')-r l4 (0)))A J 3 (^t') 3 

i 1 

Since 

ni(t')- ni (p) =-ft(0 
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we have 

9j = ^ ^'J»exp(^^^| ni= „ i(0)3l (O)A J 3 (t-t') 3 

and 

rt ^p 

9j = ^ dt'J* exp(J2 ~dt~^= n ^ 1 fe ) A K f - *') 3 

Then the last equation is the one dimensional equation with solution, which 
describes the nucleation period. This equation is the same as the in the situation 
of one component nucleation. For the first time this fact was noted in [S] , then 
it was confirmed in [S], 0, |1U| . So, there are no problems in solution and in the 
description of the nucleation period. In [9] this known fact was presented as a 
new result. We need to discuss it even in the one dimensional case because the 
errors from were accurately reproduced in |H] , which points that these errors 
were made not occasionally but were the consequence of deep misunderstanding 
of the nucleation kinetics. So, these errors have to be corrected explicitly 41 . 
The last equation can be rewritten as 

9j = w [ dt' eM-Pg 3 (t'))(t - t'f 

Jo 

with constants 

W = J* A 3 

P = - X)~^I«*=™*(°)& fc /& k > 

The physical sense of renormalizations made in |5] is wrong, the mentioned 
sense announced in |5] would be correct if all densities n.j are constant up to 
the end of the all surplus vapor depiction (due to the vapor depletion these 
densities will decrease). Actually there is no need to fulfill renormalizations 
with a physical sense, it is sufficient to write equation for Gj = Pgj(t') as 

Gj=WP [ dt'exp(-G j (t'))(t-t') m ; m = 3 
Jo 

The scale renormalization t — > {WP) X 'H leads to 

Gj = [ dt' exp(-Gj(t'))(t - t') m - m = 3 
Jo 

One can prove that this equation has the unique solution and one can get 
this universal solution numerically This approach was outlined in |15j where 
all details can be found. The base of application is the rapid decrease of J(t) 
after the end of nucleation period, which allows to present the influence of 
after the end of nucleation via the first three moments of the distribution, i.e. 

41 Solution of the one-dimensional equation in a more general case was presented in [S]. 
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via / °° t n J(t')dt' , i = 0, 1, 2, 3 (including the zero moment, which is the leading 
term). 

The iteration solution with the zero approximation Gi = and the recurrent 
procedure 

ft 

G f = / dt ' exp(-Gf (t'))(t - t') m , m = 3 
Jo 

can be also applied. But it is necessary to stress that the reason of the rapid 
convergence of iterations is the big power m = 3 in the integral term. This fact 
was noticed for the first time in [T7]. This property is the concrete realization 
of the principle of the avalanche consumption formulated in (221 - Nevertheless 
in |5] it was stated that already the first iteration gives the appropriate result 
for all positive powers m instead of 3 in the integral equation. This conclusion 
is wrong, which can be seen from the explicit solution at m — > +0, namely 

d = ln(t + 1) 

This function strongly differs from the first iteration 

g\ = t , 

which is the correct form of the first iteration for m — > +0, and also differs from 

which was used in [S]. The most surprising fact is that the main part of 9 
is devoted to get the estimates of the accuracy of the first iteration and the 
"necessary" estimate has been proven. Certainly, this estimate is wrong. We 
needn't to follow procedure [§] but can simply compare the universal 42 solution 
Gi and the first iteration G^. Both curves are J = exp(— Gi) and Ji = 

cxp(— G^) are drawn in Fig.l. and they can be hardly separated. It is very 
simple to show this approximate coincidence analytically. 

The description of the nucleation period is completed, but the period of 
consumption of the main quantity of surplus substance has to be described. 
This period is the main in the phase transition because namely during this 
period the main quantity of surplus substance is going from the mother phase 
into a new liquid phase. Really, at the end of the nucleation period the relative 
quantity of the substance is small. It can be estimated as 

where the duration to of the nucleation period can be estimated as 

42 For approximation of the universal solution see 1211 . 
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Figure 1. The form of the size spectrum and the first iteration in dimensionlcss 
units. 



(here the first channel is considered). 

This estimate is based on the analytical form of the first iteration and it is 
clear that this estimate is too long. At least the values gi satisfy inequality 

.dF e x n t 
9i < k— = 

(the critical value is referred for the first channel), which shows that the relative 
quantity of substance in droplets is small at the end of the nucleation period. 



3.1 Standard evolution after the end of nucleation 

Now we shall describe the period of consumption of practically all surplus sub- 
stance from the mother phase. Precise solution of the evolution during this pe- 
riod is unknown, the solution presented in [5] is wrong. This solution artificially 
without any comments and justifications prolongs some relations established at 
the nucleation period. This prolongation isn't correct. Here we shall present 
some effective approximate solutions of kinetic equations. Here we have no ne- 
cessity to ensure such high precision as at the nucleation period because there 
is no further nucleation (except the special case, which will be described in the 
final part). So, simply the high relative accuracy in rii is necessary. 

At first we shall use the property that at the end of nucleation the relative 
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quantity of surplus substance is small. Then it is possible to use the monodis- 
perse approximation for g^. This approximation is the following 

gi = Nvi 

where ./V is the number of droplets (It can be taken as the result of the nucleation 
period description) N = f Q J s i(t')dt' as J»io0.9. Here Vi is the coordinate of 
the monodisperse spectrum. When the monodisperce approximation for g^ is 
bad, then g^ is negligile in the substance balance. 

It is necessary now to write kinetic equations for the time evolution of v^. 
These equations can be written in approximation of the regular growth. Then 

^ - \vi t S[m(t = 0) - Nui - n t ooUS] (3) 

S = (ZGiry^V 2 / 3 , V = Y. Wi = ^ 

j 

Here and later a, are included into t . Then 
dV 1 / 3 1 

= ^(seTr) 1 / 3 t [m(t = o) - Nvi - m oc&fS] (4) 

i 

Here we conserve v\ in [rii(t = 0) — Nvi — rij oo&/<(£)]i which can lead to 
perspective approaches when for Vi some approximations are suggested 43 . One 
can rewrite the last expression as 

-jyl/3 i V -i 

— (36n) 1 / 3 J2viVit{n i (t = 0)-NZ i —=-n ioo t i fi(g)} (5) 



dt 12 



For £ the evolution equation will be the following 

H = f§?^ 367r ) 1/3 [( n «(°) - mZ - ni ooUi)v t 



3 



(6) 



It is clear that one can not not solve two last equations analytically (at least 
for the arbitrary coefficients of activity) . The analytical solution announced in 
|H] as the precise one is wrong. Now we shall suggest different approximations, 
which allow to solve these equations. 

Some useful approximations can be found in [5], [7] where it is shown that 
they can be presented as the iteration procedures with rather high rate of con- 
vergence. 

43 For example, for fast component one can put i/j to its final value (see later). 
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Approximation of ideal surface 

In J2J we shall consider S as known function of time. Suppose that S is 
known in some approximation. Then one can see the following iterative proce- 
dure 



d^ k+1) 1 



Vl t S[m(t = 0) - iV^ (fe+1) - n t oofc/id)] (7) 



dt 4 



3 

,(°) ™. *„W, „(o) 



For initial one takes Vj = A 3 i 3 . One can take in Q ^ as 

6 



or 



6 



^3 3 



(fe+i) 



The first opportunity leads to the linear differential equations, which can be 
solved and it completes the iteration step. 

Approximation of the initial concentration 

Consider (Jjjjl and put in this equation all £j to their initial values £i(0), which 
are known from the consideration of the nucleation period. Then 

jyi/3 i y 

—jT- = -:{i^) 1/3 Y, V ^ *[»*(* = °) " ^*(°)-^T - ^ «.6(0)/i(f(0))] 
eft 12 - w(?(0)) 

This equation can be easily integrated, which gives 

^ 1/3 dV^ 3 

i(367r)V3 E . ^ t[ni ( t = o) - at6(o)^ - «, oommm] ~ 1 

Approximation of the characteristic concentration 

Consider (0) and put in this equation all & to their characteristic values £i( c h) 
during the most intensive consumption of the metastable phase. The charac- 
teristic value of concentration has to be chosen as some artificial parameter, 
sometimes it is useful to choose the mean value between the initial and the final 
concentrations (the final concentration can be found without explicit solution 
of the evolution equation; it will be done later). Then 



dV 1 / 3 1 , 1/av ^ r ,. _ V 



df~ ~ ^( 367r ) t [m(t = 0) - N^ ch) -j^ — - - n» oa€i(ch) fi(€(ch))] 



This equation can be easily integrated 

r yl/3 dvv* 



la ^(367r) 1 /3^ iWiW4 t [ ni (t = 0) - {ch) Y - n, ( ch )fi(€(ch))} 
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Approximation of stationary concentration 

Equation JHJ is a typical relaxation equation, which leads to the relaxation 
of £j to the "bottom" of kinetic channel. So, it is reasonable to put & to the 
stationary value & k determined by 

= [MO)- A^-n, oofc/iH i-tiY^faW-NtiJ-nj «&/i)«t *] (8) 

i 

The stationary value depends on the value of V, i.e. fe = £j fc(V), because V 
characterizes the power of the mother phase consumption by droplets. 
Then equation for evolution of V will be the following 

dV 1 / 3 1 v1 /<.v^ , , ^ ... V 



ni fc(V0/i(£(V))] 
This equation can be easily integrated, which gives 

ctt^ 3 

o ^(36^)1/3^.^ t [ ni (t = 0) - fcOO^yy - m k(V)fi({ k (V))} 

This solution is a very good approximation to the real solution. In all sit- 
uations when there is no other big or small parameters (rather often the ratio 
fH ooVt i/nj ooV t j is a small parameter 44 , where % and j are two different com- 
ponents) this approximation works well. 

Certainly the system of algebraic equations JSJl together with ^ & = 1 has 
to be solved, which isn't a simple task. 

Generally speaking one can see that the problem is to get the true value 
for concentrations. Then the problem becomes much more simple. One can 
formulate here the following statement: 

• (F) When some approximation £j appr for the concentration is adopted 
then the kinetic equations can be integrated as 

* yl/3 dV 1 / 3 



j^(367r) 1 /3^. WiUi tn 



where 



V 

n = [tli(t = 0) - appriV) — =! - n 4 oo^i appr(V)fi(£appr(V))] 

v{£appr{V)) 

Now we shall analyze approximations, which are suitable at different periods 
of condensation. One can define initial period of condensation and the final 

44 This will be considered later. 
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period of condensation. The initial period is the period, when the quantity of 
the substance in the droplets is small in comparison with the total quantity 
of surplus substance. Initial period covers nucleation, which is appearance of 
droplets, but does not coincide with nucleation. 

The final period is the period when the quantity of remaining surplus sub- 
stance in vapor is small in comparison with the quantity of substance in droplets. 

These definitions will be slightly modified later. 

The values which correspond to the transition 45 of all 46 surplus substance 47 
into the droplets will be called the final values and will be marked by the lower 
index fin. 

Between the initial and the final periods the intermediate period takes place. 
Namely during this period the macroscopic manifestation of a phase transition 
takes place. 

Approximation of different stages 

One can analyze the characteristic features of the last solution. During some 
period (initial period) one can see that V grows rather fast but V is relatively 
small and in rj one can neglect V. Then 



1 ^(367r) 1 /3^.^ I , i t r i 



where 

h = [m(t = o)-m «,& k (v = o)m k (v = o))] 



Then 



VV3 

t 



The values of k are known - they are the coordinates of kinetic channel at 
initial values of rij, i.e. rij(i = 0). So, there arc no problems in integration of 
the last equation. The denominator doesn't depend on V. As the result V ~ t 3 . 
The next step in approximation is the following: The value of V grows rather 
rapidly. Then in there appears a rapidly growing function V, which turns 
J2i v i v i to zero 

o = v * v * * r * ■ 

i 

Then V stops to grow and in the system one can see some relaxation processes, 
which can be described as the final period. This period is characterized by some 
approximately constant and, thus, prescribed value of V = Vfi„. 

At the final period the value of V becomes the function of £ and can be 
found from the last equation, which can be written as 

= 5>iVi t [m(t - 0) - NZi(V)—¥— - n t ooUv)h{i{v))] 



45 The end of this transition is imaginary. 

46 The density of vapor equals to the density of vapor saturated over the plane surface. 
47 May be the concentration is fixed. 



30 



Later we shall write that V — V/i„(£). 
Now equations on £ can be written as 



(9) 



This equation is the ordinary relaxation equation but certainly it is nonlinear 
and multi-dimensional. The analytical solution in the general case is unknown 
but the relaxation to the stationary value, which satisfies the system of algebraic 
equations 



[MO) - Ntfi^Q - m x Ui)vt 



-Hi 5>i(°) - ^^§® ~ "3 oo^/iM s ] 



(10) 



is rather obvious. 

In the binary case there are no problems in solution of kinetic equation (JjJJ . 
It can be rewritten as equation on one concentration £i 

dfc (36tt)V3 

; — ^isiJ 



where 



= (1 - &)[(n x (0) - A^/m - n a oo6/i(a)H i - 
6(n 2 (0) - N{1 - £i)0 /in - «2 oo(l - Si)/a(l - &))«t 2] 

and in approximation i? ^„ = const this equation can be easily integrated 



(«lgl+^2(l-4l)) (367T) 1 / 3 

V 3 / 5 4 



= t 



Here appears a problem what extensive variable we have to fix in the last 
relation. When we come from V/m to i?/j n we have 



1 — #(£1) 



It seems reasonable to choose that extensive variable, which variations are small. 
The problem is what coefficients A; in extensive variable \vi we have to 
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choose (for $ all A* = 1, for V we have Xi = Certainly, we have to choose 48 
\i to have the smallest variations 49 of J^. Ai^. Here we can repeat all previous 
constructions. For the reason of physical clearness we shall conserve below V 
(and but it is more accurate to take J^. XiV%. In any case this problem is 
rather artificial because the concentrations approach to the stationary values. 

The possibility to get analytical solution in the binary case allows to suggest 
the pair approximation. We choose the fastest component (due to the kinetic 
coefficient). All other components will be the second component 50 . Then one 
can get the solution presented above. The concentration of the given fast com- 
ponent will approach to the final value. Then the same can be done with the 
set of the rest (formally slow) components. This loop can be repeated up to the 
end. 

During this relaxation such an interesting phenomena as the partial disap- 
pearance of some components from droplets can take place. Really, suppose 
that there is a fast component and a slow component. The chemical potential 
of a fast component has a minimum near £k(V = 0). The chemical potential 
of the slow component has a giant minimum far from !;k{V = 0). What will 
happen in the system? At first £ « £k(V = 0) and the density of vapor of the 
first component approaches to the value corresponding to the minimal value of 
the chemical potential for the fast component. But later the slow component 
saturates the droplets and this rises the chemical potential of the fast compo- 
nent from the minimum. Then the chemical potential of the fast component 
in the vapor has to grow. It means that the density has to grow and there is 
a partial disappearance of the fast component from droplets back to the vapor 
phase. 

One can speak about the reverse phase of the phase transition for 
some components. Certainly, this is only the first peak of the possible os- 
cillatory regime, which can take place here. This example shows us that many 
interesting phenomena were thrown out in a wrong solution presented in 
This redistibution of substances in droplets has nothing in common with redis- 
tribution of substance in droplets during the coalescense described by Lifshits 
and Slezov (3|. 

As the result one can speak about the inertial phenomena in the phase 
transition. This conclusion takes place in the multi-dimensional decay kinetics 
and is the principally new feature of the phase transition. 

We have extracted the initial period of condensation and the final period of 
condensation. We recall that the period between the initial one and the final one 
will be called the the intermediate period 51 . Do the initial and the final periods 
cover practically the whole period of condensation before the coalescense? Due 
to the avalanche consumption of surplus substance (or due to the rapid growth 
of V at the initial stage) one can establish the following very important fact 

48 With restriction ■ Ai = 1. 

49 It can be derived in some initial approximation. 

50 Or one can divide all components into two groups. 

51 For every component the positions of these periods are specific. 
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• The intermediate period between the initial period and the final period is 
rather short (in comparison with the initial period) . 

Then with a small relative error one can speak about the boundary between 
the initial period and the final period and require the coincidence of correspond- 
ing solutions at this boundary 52 . This gives the values of arbitrary constants in 
the solution at the final period. 

Now we shall discuss these periods in details. 

Initial stage 

To give the detailed description of the initial stage we shall write kinetic 
equations. It can be done as 

^ = i«, tptof'^viutffilmit = o)-m ooUid)] (n) 

i 

We see that here 

6=6 fc (t = 0) 

Then 

v X J* = U 

The same has been done in investigation of the nucleation period, but it is clear 
that this description can be prolonged on the whole initial stage (normally, 
nucleation period belongs to initial stage). One has to note that in [Sj there is 
no clear picture of the characteristic periods and stages during the condensation. 
Moreover there are many errors. For example it is stated that at the first period 
(nucleation period) the densities of components in vapor were considered to be 
constant (UJ, page 391). This conclusion is wrong. The densities are functions of 
time, the simplification is approximate equality of the densities to some constant 
values. But in consideration of nucleation, which is very sensible to the values 
of densities one can not consider these densities as to the constant values. Here 
the requirement of precision is more weak 53 - we need only the high relative 
accuracy of the densities defining. So, here we can put them to constants and 
consider the initial stage. 

One can use another possibility to enlarge the region of description by the 
analytical solution. This possibility lies in the approximation of high supersat- 
urations. Let us consider this approximation. 

In one component condensation one can show that the intensive condensation 
is possible only when 

n/n x > 1 

This conclusion comes from the standard Gibbs rule which states that the free 
energy of the critical embryo is three times less than the surface energy of the 

52 To get this boundary one has to prolong the solution at the initial period to the intermedi- 
ate period. If we prolong the asymptotic solution we can incorrectly get the long intermediate 
period. 

53 Except the situation with the secondary nucleation, which will be considered later. 
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same embryo. If n/n^ — ► 1 then _F C is too big and the nucleation rate, which is 
proportional to exp(— F c ) is negligibly small 5 . Ordinary n/rjoo w 3 -j- 5. 

In the situation of multicomponent condensation this requirement {n/n^ 3> 
1) isn't necessary because there are many components. But still there can 
exist some leading components, for which rij/rij ^ 3> 1. Then one can neglect 
n i oo£,i.fi (£,) in comparison with m — Nvi and then get 

^ = t S[m{t = 0) - Nvi] 
Analogously one can get a more simple equation on & 

^ = ^ll( 367r ) 1/3 [(^(0) - i ~ 6 £(%(0) - (12) 

j 

Then the equation on the coordinates of kinetic channel will be 

[MO) - mZ) Vt i 6^(^(0) - a&^h = o 

Then the value 

(m(0) -N&-)v t iHi = c 
v 

is some invariant (it doesn't depend on component). 

It is clear that the channel will be purely kinetic. It is easy to see the 
concentrations & can not approach to 0, because then the equality 

V 

ni(0)=^(c/v ti + N-) 

will be violated 55 . 

One can easily get equation on c. Really, 

n»(0) 



6 



c/v t i + N\ 
Then one can get equation for v via c 

m(0) 



'c/v ti + N^ 

Then v — v(c) and because of ^ = 1 one can get equation on c 



1 



^c/vti + N^ 



54 The value of surface tension referred to one molecule can not be too small otherwise one 
has to take into account the fluctuations of the surface and the situation will be close to the 
second order phase transition. So, here there is a hidden contradiction between the intensive 
nucleation and the possibility of thermodynamic description of the critical embryo. 

55 When v is the regular function of concentrations. 
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This equation can be easily solved. Ordinary it has a unique solution. 

As the result we see that to observe the kinetic channel near £j = it is 
necessary to have really giant values of It means in this case that the values 
of £i are governed by fi and they are rather insensitive to the vapor depletion. 
Then & k is the quasiintcgral of evolution and the stationary approximation for 
& k suggested earlier is good. Then it is easy to see the relaxation to this value. 

Final stage 

The final stage can be described by approaches described earlier. Here we 
shall see that the end of this stage can be described extremely simple. The 
final values (they can not be attained because of coelescense) can be gotten as 
solutions of the algebraic system 

di>i/dt = 

or 

o = m{t = o) - Nvi - m oo6/i(l) (13) 

They will be marked as Vi fi n . When z/, are close to Vi f in , then the kinetic 
equations for At^ = Vi — Vi f in can be rewritten as 



dAv; 1 



dt 4 

where 



-V t iNSfinAl/i 



S fin = (36n)y^v^ k) f J 



2/3 



These equations are linear and can be easily solved. As the result we have the 
solution at the end of the final stage. The solution is 

a , VtiSfi n Nt 

l\Vi = Cj exp( ) 

with generally unknown negative constants Cj. They can be determined by co- 
incidence of solutions at initial stage and at final stage at the boundary between 
these stages. 

This completes the general investigation of the one-stage nucleation. 

In the next section when we shall study the multi-stage nucleation process. 
We shall give another approximate solutions, which will be valid under some 
specific situations of small concentrations. 



4 Nucleation in the case of small concentration 
in droplets 

Why shall we pay a special attention to the further investigation of nucleation? 
One can recall one important feature of one component nucleation. The one 
component theory predicts that the number of droplets N is proportional to the 
amplitude of spectrum / in the power 3/4. Since / is proportional to exp(— AF C ) 
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it means that for two channels with AiF c < A 2 F C — 1 the strong inequalities 
/l ^> h and N\ » N 2 take place. This conclusion in the multicomponent 
condensation isn't valid because one has to take into account specific kinetic 
factors which can be different for different channels. 

One can also notice another important feature. In the multicomponent the- 
ory the positions of thermodynamic channels don't observe the order. When in 
the one component case we analyze the dependence F($) along the thermody- 
namic channel we sec that the higher is activation barrier the higher lies the 
curve (along the bottom of thermodynamic channel). But since the values 
of surface tensions in different thermodynamic channels are different one can 
come to situation when the channel with the highest barrier (i.e. wiht the high- 
est F c ) later becomes the lowest channel 56 . It means that we have to analyze 
all channels despite the heights of their barriers. 

At first we shall describe the ordinary situation. The process of nuclcation 
is ordinary localized in time. Really, already in the initial period 

Vi = A 3 t 3 

and 

m(t) =n i {t = 0)-NA 3 i t 3 . 

Thus, every component leads to the fall of intensity of the droplets formation 
in the channel number j (let us enumerate the channels according the heights 
of activation barriers) as 

J{t) = J(t = 0)exp(^% c i jin^-mit = 0))) = J(t = 0)cx P (-^z/ c i jNA 3 t 3 ) , 

i i 

where v c j j is the number of molecules of i component in the critical embryo for 
the channel number j. This rapid decrease of the nuclcation rate J according to 
the previous relation means that there will be no further formation of droplets. 
So, the process of intensive nucleation is finished in all channels. This situation 
has been already described in the first part of the previous section. 

But still there is one important exclusion from this situation. We shall con- 
sider the nucleation in some channel which will be called as the second one. So, 
the last index now will merk the channel of nucleation. When the concentration 
£i of i-th component in the droplets in the first channel is small, then Aj can 
be small and at moderate times t the value — v c j i{n(t) — m(t = 0)) is small 
enough when v c i 2 isn't too big 57 . So, the intensive nucleation (in compari- 
son with initial value) in the second channel continues even when the intensive 
nucleation in the first channel stops. Then the number of droplets appeared 
in the second channel becomes important and it is interesting to calculate this 
number. 

The reason of termination of intensive nucleation in the second channel 
can be the vapor depletion by both droplets appeared in the first channel and 

56 We compare channels 1? being fixed. 

57 We suppose this property for the second channel. 
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droplets appeared in the second channel. This situation becomes complex and 
requires a special investigation. It will be given in this section. 

From the first point of view it seems that this situation is rather artificial. 
But as has been already noticed in introduction when we have two immiscible 
substances we come to two separate processes of nucleation - the first process is 
the nucleation of the first substance into the droplets of the first substance and 
the second process is the nucleation of the second substance into the droplets of 
the second substance. These processes have nothing in common. This situation 
has to be described in frames of the binary condensation also and namely this 
situation falls into the case, which will be considered. When there is a small 
miscibility of substances we have a moderate influence of the first process on 
the second one and have to take this interaction into account. 

We shall call the situation when several channels of nucleation are sequen- 
tially important as the "long nucleation". 

Now we shall investigate the case of long nucleation in details. We shall 
describe here the two component case, then the necessary remarks for general- 
ization will be given. 

We have already noticed the possibility of partial reverse condensation and 
shall call this as "inertial" effects. 

At first we shall analyze the general features of the process, then we shall 
give the theory for the process with no "inertial" effects in the first channel, 
then the situation with "inertial" effects in the first channel will be described. 



4.1 General features of the "long nucleation" 

The nucleation process in the first channel has been already described and we 
know the number N± of the droplets appeared during the nucleation in this 
channel (earlier this quantity was denoted by N). Let us suppose that in the 
first kinetic channel £1 k is small. In the second channel £2 c is small. One can 
justify the following statement: 

• Up to the end of nucleation in the second channel the droplets appeared 
in the second channel can not essentially perturb the rate of growth of 
droplets appeared in the first channel. 

Then we can investigate the separate process of vapor depletion by the 
droplets appeared in the first channel. 

For the evolution of droplets appeared through the first channel £1 <C 1, 
£2 = 1, the values of concentrations arc known. 

The system of kinetic equations for this problem can be written in approxi- 
mation of ideal solution since the solution is really the dilute one. The system 
of equations is given by 58 

^ = ^S(ni(0) - ATi/i -m oofr) , 

58 For dilute solution one can use the approximation of the ideal solution (may be with 
renormalized parameters). 
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— ^-S(n 2 (0) - Nv 2 ~ n 2 oo) , 

S = (36Tr) 1 / 3 (v 1 v 1 +v 2 V2) 2/3 ■ 

The value of S when there is no extremely small values in ratios Vi/vj can be 
presented as 

S=(Mir)^(v 2 v 2 ) 2 / 3 . 
Then it is possible to integrate these equations 59 . Equation for v 2 will be 

±1 = ^(36^/3(^)2/3(^(0) - Nv 2 -n 2oo ) . 

Since the r.h.s. doesn't contain t this equation can be integrated 

P dv2 

So, now v 2 is the known function of t. Then equation for v\ can be written as 

^ = ^i(367T) 1 / 3 ( U2 ^ 2 (i)) 2/3 (n 1 (0) - Nv x - ni oofi) 
and since £1 w v\/v 2 it can be rewritten as 

^1 = ^J:(3 67r )l/3( U2 ^ (t)) 2/3 (ni(0) _ N _ ^/^(t)) . 

at 4 

Since it is the first order linear differential equation it can be easily integrated 
by the known formulas. 

This solution is explicit, but the model with the small concentration and 
ideal solution can be criticized. One can see that the limit of the last solution 
is 

m(t) -► . 

The final values of u n are obtained on the base of algebraic equations l|13fl 
and they doesn't depend on v t t and on ratios 60 of n, oo. So, one can easily 
come to the situation when £i /j„ doesn't go to 0. The last solution doesn't 
satisfy this condition. 

The last remark means that we have to investigate the general situation. 
From the previous section we know that the solution can be composed from the 
part corresponding to the initial period and the part corresponding to the relax- 
ation during the final period. So, there appeared two characteristic problems: 

• 1). To describe nucleation in the second channel during the initial period 

59 In the first equation one can neglect ni oo£i, which leads to 

at 4 

and the first equation can be then integrated in a very simple way. 

60 The initial positions of kinetic channels strongly depend on vt i and ri; oq. 
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• 2). To describe nucleation in the second channel during the final period. 
The last problem has to be solved in different situations. When 61 

£l fc/l(£l fe) > £l finfl(€l fin) 

then ordinary there are no inertial effects 62 . The rhs in the last inequality can 
be taken at initial values of vapor densities since as it has been noticed the 
channel in the situation of small concentration is very deep and concentration is 
here a quasiintegral of evolution. One can see here that the value n\ constantly 
decreases in time. Here one has to investigate the nucleation in the second 
channel during the relaxation with decrease of density. 
When 

£l fc/l(£l fe) < £l finfl(€l fin) 

one can come to the situation when m decreases at first in time to the value 
£1 fc/i(£i fe)- The value of m can not come too close to this value. Then it 
begins to increase and approaches to £1 /m/i(£i fin)- So, beside the problem to 
describe the nucleation during the decrease of m it is necessary to describe the 
nucleation during the increase of n\ . This problem also splits into two separate 
problems to describe the evolution when there is essential increase of n\ and to 
describe evolution when m already attains the asymptotic value. 

It is extremely important to note that we need the solution only when £1 <C 1. 
Only when d«l we can suppose that there will be rather intensive droplets 
formation in the second channel. At the same time only then the solution will 
have a very high precision, which is necessary to describe the nucleation in the 
second channel. 

4.2 Secondary nucleation during the initial stage 

We shall call the appearance of droplets through the second channel as the 
secondary nucleation. During the initial stage of the primary condensation the 
equation for the density n\ as function of t will be 

ni(t) = ni(0)-JVi(Ait) 3 -<& x . 

Here N\ is the total number of droplets appeared through the first channel (since 
the nucleation through the first channel is over we know this number), g 2 i is the 
number of molecules of the first component in the droplets appeared through 
the second channel. Then having repeated for the second channel the same 
procedure as we had done for the first channel we get the following equation for 
the number of the molecules of the first component gi i in the droplets appeared 
through the second channel 

92 i(t) = .h *A? 2 f dt'(t - t') 3 cxp(-[^ „. 2 c ( T )N 1 A*t' 3 + 

Jo_ 

61 The final values here and later are taken only with account of droplets appeared through 
the first channel (this number is already known). 

62 Certainly, the opposite situation can take place. 
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y^Vj 2 c(r)£i k92 1 (*')/&■ k]) 

i 

where t = or r is some characteristic value, may be at the asymptotic solution. 
Here Ai 2 is the constant in the law of growth 

1 1/3 



of the molecules of the first component in the droplets appeared through the 

second channel, J 2 * is the amplitude value of the stationary rate of nucleation 

in the second channel, ^ 2 c is the number of molecules of i-th component in 

the critical embryo in second channel, Aj is the constant in the law of growth 
1/3 

dv i 1 /dt = Ait of the molecules of the i-th component in the droplets appeared 
through the first channel, & k is the concentration of i-th component molecules 
in the second thermodynamic channel. 

Having combined the constants we have 

92 i(t) = Ji\ dt'it - t'f cxpi-At' 3 - Bg 2 lit')) 
Jo 

with known constants J 2l A, B. Both functions exp(— At' 3 ) and exp(— Bg 2 lit')) 
can be interpreted as step-like functions, which allows to state that 

To get the rate of nucleation the value of g 2 1 in the r.h.s. can be calculated 
approximately without account of exp(— At' 3 ) in frames of the first approxima- 
tion in the standard iteration procedure 

92 1 (0) =0, 52 1 (i +i)(t) = J2 I dt'it - t'f exp(-A*' 3 - Bg 2 x (i) it')) 

Jo 

So, <? 2 1 (1) = J2^ 4 /4 and then 

52 1 (2) = J2 [ dt'it - t'f cM-At' 3 - BJ 2 t' 4 /A) 
Jo 

So, the rate of nucleation looks like 

J = J 2 exp(-^' 3 - BJ 2 t' 4 /4) 

and it has the step like behavior. Moreover it is the product of two step-like 
functions exp(— At' 3 ) and exp(— B J 2 i /4 /4). So, we can introduce two character- 
istic sizes 

A a t = A- 1 ' 3 

and 

A b t = (BJ2/4)- 1 / 4 

We have to take the minimal value from A a t and A^t to determine the nucleation 
duration 

At = min{A a t, A b t} 
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Then the rate of nucleation as a function of time is 



j(t) = j 2 e(t)e(Ai - 1) 



and the total number of droplets N 2 appeared through the second channel is 



Now we have determined all main characteristics of phase transition. These 
characteristics can be determined with the more high accuracy using the per- 
turbation technique and some iteration schemes. Here it is important that the 
base for these constructions is now given. 

5 The properties of further evolution and their 
effects on the secondary nucleation 

Now we have to consider the situation when the secondary nucleation takes place 
at times greater than the initial period. We know that at some intermediate 
times after the end of initial period but before the asymptotic final period the 
solution is known with essential error. Fortunately one can prove the following 
statement, which will help to give a simple analytical description: 

• When the process of the secondary nucleation stops after the end of the 
primary nucleation initial period but before the time, which two or three 
times exceeds the time of the end of initial period of the primary nucleation 
one can simply neglect the vapor consumption by the droplets from the 
second channel after the end of primary initial period. 

Here we define the initial and the final periods referred to the correspond- 
ing components. 

Now we shall explain this statement. At first we shall see that the period 
between the initial period and the asymptotic period has a rather small (in 
comparison with initial period) duration. We shall see this on example of the 
one component nucleation where the precise solution can be easily obtained. 
Really, the evolution equation can be reduced after the certain renormalizations 
to the following form 

dz 3 
— = 1 - z A 
dt 

with the initial condition z(t = 0) = 0. Here z is the renormalized coordinate 
of the droplets linear size -i? 1 / 3 . The solution is the following 





and the integral can be taken analytically, which gives 




1 1 
— arctan — — 
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Figure 2. Different approximations to the evolution of density 



The last solution has the asymptotics 
1 1 — z 

t = — In — = = z = 1 — exp(v37r/6)v3exp(— 3t) 

3 V3cxp(V3^/6) V 1 ' V ' 

at big t and 

t 4 3 i 7 
4 4 7 

at small t. 

In Fig. 2 one can see the precise evolution and the asymptotics. The curves 
are drawn for the value of excess of density, i.e. for the normalized difference 
^ ~ n — rioo, which is proportional here to 1 — z 3 . One can see one long curve, 
which binds the curves, which are cut off. This long curve is the precise solution. 
The short curves are asymptotics. There is one asymptotic f as for the big times 
and two asymptotics for the small times. The asymptotics for the short times 
differ by the number of terms for z in l|14f> . which are taken into account. The 
curve /m corresponds to the account of the first two terms in expression for z 
in HI 4(1 and the curve /( 2 ) corresponds to the account of the first three terms in 
expression for z in 1 (14(1 . 

One can see that the period of time when both the long time asymptotics 
and the short time asymptotics don't work has a rather short duration. This 
duration has to be compared with the duration of initial period. Since the 
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precise solution is known one can give the estimates but everything is clear 
already from the picture. One can only note that this property comes from the 
big power 3 in the law of droplets growth. This is the same reason as for the 
property of the avalanche consumption |22| . So, one can here note the following 
important property 

• The period of consumption of the main surplus substance is rather short 
in comparison with the initial period. 

So, one can speak also about the avalanche consumption of the all surplus 
substance. 

This property is typical also for the case of the dynamic conditions |17|. 
Earlier it was only noted that the period of consumption of all surplus substance 
is short in comparison with the period of waiting (see |17|). Now we can state 
that the period of the avalanche consumption (or the intermediate period) is 
short in comparison with the period when the vapor density approximately 
equals to the density at the moment of the most intensive formation of droplets 
(or simply to the density at the moment of intensive formation of droplets), i.e. 
to the initial period. 

The small parameter characteristic for the initial period in the solution of 
the last equation is g = z 3 in comparison with 1. One can see that when 
771 3> 1 the value z m will be small even when z is comparable with Zfi n = 1 
(for example (1/2) = 1/8 <C 1 can be considered as the small parameter). The 
small parameter for the asymptotic decompositions is 63 

|z-Z/«n| < l 
Zfin 

where zji n corresponds to Vfi n . Here Zfi n — 1 and it means that simply z is 
close to 1. So, we see that there is practically no gap between the initial and 
the final asymptotic periods. 

One has to outline that there is no satisfactory approximation for the value of 
the vapor density or for the value of the number of molecules in the liquid phase 
as a function of time because namely during the intermediate period occurs the 
phase transition - the main quantity of substance goes from the vapor to the 
liquid phase. Fortunately, these values aren't directly involved in the secondary 
nucleation description. But for the times of duration of characteristic periods, 
which are necessary for the investigation of the secondary condensation the 
required approximations are rather precise. 

In the multicomponent case the convincing arguments are the same. 

Now we shall explain why it is possible to miss already the two or three times 
of duration of initial period. Certainly the density attained at the asymptotics 
n t as is essentially smaller than the initial density. So, the rate of nucleation 
will be many times smaller than the initial rate of nucleation. The value of the 
critical number of molecules of given component in the second channel can not 

63 This estimate is very rough, it depends on activity coefficients /; and can be established 
only in the order of magnitude. 
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be very big (since the derivative of the free energy barrier on density is this 
value) otherwise the rate of nucleation at rij as will be absolutely negligible in 
comparison with the initial rate of nucleation. Then the period of initial peak of 
nucleation through the second channel will be not negligibly short in comparison 
with the duration of initial stage. Then we see that the number of droplets 
though the second channel appeared at the asymptotic can be comparable with 
the number of droplets through the second channel appeared in the peak of 
nucleation only in one case - when the duration of the whole nucleation period 
(including the nucleation at the asymptotics) strongly exceeds the duration of 
the initial period. 

A special question which has been dropped out of consideration is the mul- 
tiplicity of the final states. It was supposed that the "final" is unique. But the 
algebraic equations appeared from the conditions of the absence of the regular 
growth can have several roots. It means the several possible final states 64 . This 
complicates the theoretical description. The first problem to solve here is to 
see what final state corresponds to the given channel of nucleation. We have to 
construct asymptotic for every final state and then to decide what asymptotic 
is the closest one to the solution at the initial period, this problem is the pure 
algebraic one and can be solved without difficulties. 

The second problem is the transition between final states. Since the regular 
growth is absent here, we can return to the free energies of the final states (or of 
the droplets in the final states) Ffi n . The rate of transition between final states 
is proportional to exp(— AFfi n j j) where AFjj„ j j is the height of transition 
barrier (it depends on the numbers i and j of final states). The values AFfi n i j 
are so giant that ordinary there is no sense to consider transitions between final 
states. But in some peculiar situations when the height of transition barrier 
is't too big, one can observe a very interesting kinetic behavior, which forms 
the problem for a separate publication. This transition may influence also on 
kinetics of coalescence 65 , which is also a very interesting point of investigation. 
One can only note that transition between final states situated on opposite sides 
of metastability gap has to go through the total evaporation of embryos. This is 
very slow process, its characteristic times strongly exceeds even the character- 
istic coalescence times. Generally speaking, kinetics of evolution of final states 
resembles in mathematical structure of equations the micellar cvoltuion and, 
thus, we need not to give here the special description. 

5.1 Secondary nucleation during the relaxation period 

We need now to investigate the evolution at big times, corresponding to the 
asymptotic solution. 

Here we have to take an explicit account of the initial peak of metastability 
by introduction of total number of droplets appeared through the second channel 

64 " Final" states aren't the real final ones, the process of coalescence will follow the process 
of the regular growth of monodispcrsc spectrum. 

65 Certainly, the coalescence of ground final state (with minimal Ff in ) will be described by 
the classical Lifshic-Slezov asymptotic. 
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during this peak. This number is given by 

poo 

A* = / J 2 exp(-At' 3 )dt' = 0.9 J 2 A~ 1/3 , 
Jo 

where this approximation and all parameters are described in the previous sec- 
tions. 

We can here also introduce the duration of nucleation peak At as 

At = A- 1 / 3 = A a t . 

At these big times one can write for g 2 1 the following expression 

fl2 ! = N*z 3 + f 2 (z-Sz) 4 /4 . 

Here f 2 — J 2 s ({^j /m})/Ai 2 ({rii /m}) is the stationary rate of nucleation 
J% s{{ni fin}) through the second channel attained at the asymptotic values 
(it can be calculated by known formulas for the stationary rate of nucleation at 
known final asymptotic values) divided by the rate of droplets growth Ai 2 ({rii /m}) ; 
the value 5z is the shortage due to the fact that initial period is already ac- 

1/3 

counted (this value can be also explicitly calculated, it is coordinate v^' at- 
tained by the droplets born at t = up to the end of initial stage), z is the 

1/3 

coordinate v l 2 attained by the droplets born in the second channel at t = up 
to the current time t, i.e. z = f Q Ai 2 {{rii(t)})dt — > Ai 2 ({rii fi n })t. 

Since at asymptotic the concentration is constant and the densities are also 

1/3 

constant, then the values v i 2 grow in time with constant velocity independent 
on Vi 2 . Here Vi 2 is the number of molecules of the i-th component in the 
droplets born in the second channel. This explains the powers 4 in the previous 
expressions. Sometimes it is worth to add the shift 5'z to the coordinate z in 
the term N*z 3 due to the irregularity of droplets velocity at the initial stage 
and at the transition from initial stage to the asymptotic solution (when z is 
defined as Ai 2 ({"i fin})t). 

Since z = A x 2 ({m f in })t and 5'z = / (Ai 2({n;(t)}) - A x 2 ({^ fin}))dt -> 
8' z — J °°(Ai 2 ({rii(t)}) — Ai 2 ({rii fi n }))dt, the expression 66 for g 2 i will be the 
following 

92 i = + S'z) 3 + f 2 (z + S'z - Sz) 4 /4 . 

We see that the expression for g 2 j is a very sharp function of z. One can 
then speak here also about the avalanche consumption of the metastable phase, 
which was outlined in |22| as the main feature of nucleation. This property 
allows to speak about the upper boundary. To introduce the definition of this 
boundary it is necessary to recall the form of spectrum as function of time 

f = h 6xp(- ^2 Vi c 2 finii fin92 l/£l fin) ■ 
i 

Again the avalanche character of vapor consumption can be proven. 

66 The value of 8z has to be reconsidered in the same manner. 
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Then the boundary AfZ is defined according to 

Vi c 2 fin£,i finff2 l/£l fin = 1 

i 

or 

(J2 V * fi»& = ( N *( Z + <5 ' Z ) 3 + ^ Z + 5 ' Z - Sz f)/Zl f*n ■ 

i 

This equation gives the value of 2, which will be the boundary AfZ of spectrum. 
The spectrum will be approximately 

f = f 2 Q(z)Q(A f z - z) 

and the number of droplets is 

N 2 = iV* + 0.9/2A/Z w iV* + / 2 A/z . 

Ordinary the values <5z and <5'z can be neglected here. 

This completes the general investigation of the secondary nucleation but 
below the more accurate methods will be presented. 

5.2 Description of the secondary nucleation at the de- 
creasing asymptotic 

During the asymptotic (final) period the concentrations & of the droplets ap- 
peared in the first thermodynamic channel will approach the final values ac- 
cording to the relaxation equation. Namely, the final values ^ satisfy the 
system of equations 



[(rij(O) - N ^f™ _ n . ^g. finfafein^vt 



(15) 



6 fin _ N & fin Vfm ^ fl ^ ~ Tlj oot-j finfj(£fin))Vt j] = 

3 

Then equation Q can be rewritten near = £j as 

^ = iiK/m)^ - & /in) (16) 



where 



and 



L _5S 



s = ^l(36,nnm-mi^ 



(17) 



o€ifi)v t i - &^2(nj(0) - N£j Vfm _ ^ - rij ooCj/iM j 
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Solution of the last system of equations is known 
& - & fin = C l exp(- — — ) , 

<5i r t = /i„)| 

with 67 constants Cj, which can be found from the precise solution of © by 
transformation to asymptotics. These solutions lead to analogous exponential 
solutions for n, — fj n , namely 

rii - rii f in = C n i exp(- - - ) 

0% n rt 

with constants C n i and corresponding relaxation times 6i r t. 

Now we know rather precise asymptotics and can write the equation for the 
nucleation rate in the known exponential approximation 

J(t) = J2 s (7ym)exp(y~V t c fin [C„ iexp(- - ) - & fin 92 3 ]) . 

Here the lower index in c ji„ means that the critical value of Vi is taken at 
the final values of rii. The value of gi i is the number of molecules of the i- 
th component in a liquid phase appeared through the second thermodynamic 
channel. 

At first we shall analyze the case of decreasing asymptotic when even without 
the natural decrease (due to gi j) of the rate of nucleation this value decreases 
in time. 

Now we shall hold only the leading 68 exponent, which brings the last equa- 
tion to 

J(t) = J2 /i„exp(^o c fin\C n |exp(- - - )) 

Oq n r l 

CXp( ^ ^ U% c fin^i fin~Z ' ) 
i Si fin 

J2 fin = J2 s(nfi n ) , 

where we have marked the leading exponent by index (sometimes there are at 
least two leading exponents with the same times of relaxation since the sum of 
concentrations equals to one, so the amplitudes have to be summarized). Then 

J{t) = J2 /mexp(i/ c f in \C n |exp(- t ))cxp(-A,-ff 2 j) 

00 n rt 

where the constant Aj is given by 

Aj — ^ ^ Vi c 2 fin^i fin~Z 

i fin 

67 If the relaxation really takes place, all L; are negative. 
68 The exponent with the greatest relaxation time. 
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Since g 2 fi n is invariant (it doesn't depend on j), the rhs of expression for J 
doesn't depend on j also. This approximation can be used for further analysis. 

1 /3 

At the asymptotics one can write v i ' 2 = A; 2 (t — t') where t' is the time of 
formation of the droplet, 

A, 2 = (l/12)(367r) 1/3 u 2/3 4~^ / ^ m w 4 t (m as - m ^ k f in fi{^k fin)) , 

£i k fin is the equilibrium concentration corresponding to the densities at the 
asymptotic for the droplets appeared through the second channel, rn as is the 
density at the asymptotic. 

Now one can write equation for g 2 i as 

f z t' 
92 i = N*(z+5'z) 3 +f 2 j / dx(z- x) 3 cxpOo c i fin\C n o|exp(-- -)) exp(-A J g 2 j) 

JSz^a °0 n rt 

where we specify that the nucleation is going through the second channel. Here 

1/3 

z = Ai 2 t, x = z - i// ,t' = x/Ai 2 , f 2 i = J 2 S {nfin)/Ai 2 . Now this equation 
will be solved. It is easy to note that 

r % t' 

/ dx(z - x) cxp(z/ c 2 fin\C n o|exp(-- -))cxp(-A i 32 j) 

JSz "0 n rt 

grows very rapidly (faster than z 3 ). Then the spectrum of sizes is 69 

t' 

f = h j<d{z)e(S z - x)cxp(v c 2 fi n \C n [exp(- )) 

plus the monodisperse peak from the initial period and z(t) = Ai 2(^ — ^0) where 
to is the time shift (it has to be associated with C„ analogously to 5z, S'z). 
Here 5oz is the root of equation 

gj 2(S z) = l/Aj 

(j isn't important here). The total number of droplets is approximated by 

f Saz x 
N 2 = A* + f 2 j / dxexp(^ c 2 fin\C 0n \ exp(--r — -z -)) 

JSz A i 2^0 n rt 

or after linearization 

S z S z 

N 2 = N*+(5 z-5z)f 2 j+f 2 ji/ c 2 fin\C n \5 „ r tA t 2 [-exp(-- — -)+exp(- 



Aj 2 5q n r t Ai 2<5q n rt 



Here to solve equation on Sqz it is necessary to calculate g 2 (z) in appropriate 
form. One can done the following transformations: at first we shall present g 2 

as 

f z t' 
92 i = N*(z+5'z) 3 +f 2 j / dx(z-x) 3 (-l+exp(^ c 2 fi n \C n | exp(-— — -) exp(-A i .g 2 j)) 

JSz "0 n rt 

69 Ordinary it is possible to neglect Sz, but even with <5z taken into account the derivation 
will be quite similar. 
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+/ 2 ,(z-<5z) 4 /4 
Then we can out the upper limit in the integral to oo and get 

92 i = N*(z+5'z) 3 +f 2 i / dx(z-x) 3 (-l+cxp(v a c 2 fin \C n |exp(- — - — -)exp(-A,g 2 j)) 

JSz "0 n ft 

+f 23 (z~Sz) 4 /4 

Then 172 i is the polinomial of the forth power and we know the analytical 
structure of 92 i- Then the value of SqZ is solution of the algebraic equation of 
the forth power. It can be easily solved. 

When we take into account only the leading term / 2 jZ A /4 we come to the 
method described earlier. 



5.3 Secondary nucleation at the increasing asymptotic 

The rate of nucleation imaginary calculated without droplets appeared through 
the second channel will be called the ideal rate of nucleation. 

All methods developed in two previous subsections can be used here. But 
in the case of increasing ideal rate of the droplets formation through the second 
channel there is a following danger - the increasing density rii plays the role 
analogous to the increasing ideal supersaturation in case of dynamic conditions 
[T5) . So, at some density it is possible to see the compensation of increase of 
density due to evolution of droplets appeared through the first thermodynamic 
channel by the molecules consumption by the droplets appeared through the 
second thermodynamic channel. So, here we can see the peak of intensive for- 
mation of droplets through the second channel much earlier than we attain the 
asymptotic. We has to investigate this situation. 

The next difficulty is connected with the fact that due to increasing density 
the term N*(z + S) 3 can not be now considered as the reason for the practically 
instantaneous cut-off of the spectrum. This radically completes the considera- 
tion. 

We can introduce the natural length of spectrum A nat z corresponding to 
the given vapor densities rii fi„. This length is the length of spectrum when 
the initial vapor densities are equal to rii /„. It is shown in |TB] that the 
characteristic time of the spectrum formation in the situation of decay and in 
the situation of dynamic conditions will be approximately one and the same 
(see the law of inertia of characteristics of the phase transition, page 92). So, 
we have here some natural characteristic of duration of the spectrum formation. 

The characteristic time A nat t will correspond to this length of spectrum. 
One can compare A nat t with the relaxation length A re it = So n rt going from 
the asymptotics. 

When A nat t » A re it one can easily use the approach of the previous 
sections, there is nothing to change in the previous approach. 

When A nat t << A re ;t one has the situation described as dynamic conditions 
with the source, which can be effectively linearized (see ^H|)- So, here the 
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description is also known (the only specific thing is that we have to add the 
quantity N* to the number of droplets obtained at the second peak of densities 
(the first peak is the initial peak)). 

When A nat t ~ A re it one has use the modified Gaussian method, which has 
been already used to investigate the nucleation on several types of heterogeneous 
centers in dynamic conditions (see [TU]1. 

The method similar to the mentioned method will be described below. We 
shall present the simplest version which requires the minimum formulas and 
calculations (but isn't too precise). We shall describe here the nucleation in dy- 
namic condition with essentially non- linear behavior of external supersaturation 
(effective external conditions) in time (for external supersaturation see 

We shall start with the case when the ideal supersaturation can be linearized. 
The integral equation, which has to be solved can be written in the following 
form 

g(z)=a (z - x) 3 exp(c($(x) - $*) - bg(x)) . 



Here a, b, c - are some constants, $ can be interpreted as rij, index * marks the 
value in some characteristic moment t„, which is the moment of maximum of 
c($(z) — bg) and z(t*) — 0. The values z and x have the same sense as in the 
previous section 70 . 

When $(z) — $* can be linearized the equation after the certain renormal- 
ization can be written 

g(z) = a! I {z — x) 3 exp(c'a; — b'g(x j) 



with new values of parameters c',b'. 

The last equation can be solved by the following approximate procedure. 
We know that the iterations defined according to 

ff(o) = , 

g{i+x){z) =a! I (z- a;) 3 exp(c'a; - b'g (€) {x)) 

J —oo 

converge rather fast. Already the spectrum with g as the first iteration ensures 
the correct number of appeared droplets 71 (a relative error is less 72 than 0.15). 
Why it takes place? Because the droplets appeared at the first moments of 
time (when bg was really small) are the main consumers of vapor during the 
nucleation period. This takes place due to the big power 3 in the integral term. 
So, there exists a boundary z = —bo, which have the following properties 

• (') The droplets appeared before z = —bo govern the evolution during the 
nucleation period. 



70 The value <t>* can be taken as the asymptotic (without droplets appeared through the 
second channel) value (without any correspondence to the time 1211 ). 

71 Here it is taken into account that d 2 Q/dt 2 < 0. 

72 One has to see this error when $ — <!>* can be linearized and then to prove that for 
d 2 <$>/dt 2 < the error will be smaller. 
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Figure 3. Precision of the method in the linear case 



• (") The inequality bg(-bo) <C 1 takes place. 

It is clear that bo takes the value, which has the order of c' _1 (more accurate 
it is (0.7^0.8)c" 1 ). 

So, we can suggest the following procedure: At first to calculate the value 
of g defined as 

/-bo 
dx(z - x) 3 exp(c($(z) - . 
-oo 

Then to calculate the number of droplets as 

/ — b ^oo 
dzexp(c($(z) -$,)) + a / dzexp(c($(z) - - bg(z)) (18) 
-oo J — bo 

This method is rather precise, the ratio N appr /N where N appr is the number 
of droplets calculated in l|18|) and N is precise solution is drawn in 73 Fig. 3 
Now we shall apply this method to the case under consideration. 
The equation we have to study is 

92 l = h l / (z - x) 3 exp(-is c 2 fin\C \exp(-x/Soz))exp(-Aig 2 i)dx , 

J —OO 

73 Here the linear case is considered. 
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where 8$z = So n 2^Ai 2 and /2 1 has the sense J2 fi„ divided by the rate of 
growth Ai 2 - Here C ~ C n . So, here $ - $* ~ -z/ 2 c /m|Co| exp(-a;/(5 z) 
and, thus, 

d^/dx > , d 2 <5>/dx 2 < . 

One has to add the term exp(— vj c 2 ^ * — N^z* + x) 3 ) due to the 
initial peak of droplets. Here z* is is the coordinate of this monodisperse peak 
at t = t*. Then 

52 1 = /2 1 / (-2 - a;) 3 exp(-^ c 2 /in I Co I exp(-z/<5o-2)) 
</ — 00 

exp(- i/j c 2 ^ ^ — N*(z* + x) 3 )cxp(-A 1 g 2 i)dx 

Sl A: fm 

and inequality d 2 <&/dz 2 < also takes place. 
Now 

c(*(a;) - $) = -^0 c 2 /m|C | cxp(-x/M -^fjc 2 ^ - f — N*(z* + xf 

. Sl k fin 

Due to the last inequality the method will work even better than in the 
linear case because here the properties (') and (") will be better. The problem 
is whether one can calculate g analytically. 

Actually g is the polynomial in z of the third power 

s = XX(-i) 3 -v 



i=0 
-b 



i!(3-i)! 



/-oo 
x (3_i) exp(c($(z) - $»)) . 
-00 

Now we shall calculate TV and have to take definite integrals. It can be done 
by the quickest descent method since in the subintegral expressions there are 
only exponents. 



6 Discussion 

The theory for multicomponent condensation can not be presented in such a 
concrete and compact form as in the case of one component condensation. This 
is a result of the presence of complex and generally arbitrary coefficients of 
activity /j. Here we have supposed only some general properties of /, such as 
the regularity of solution at small concentrations 

The condensation in the binary systems is much more interesting than in 
the one component case. The properties of reverse condensation, of the sec- 
ondary nucleation are simply absent in the one component case. Now they are 
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predicted only theoretically and it will be very interesting to see them experi- 
mentally. Unfortunately, even the experimental results used in [5] don't belong 
to the kinetics of the whole phase transition, but to the stationary rate of nu- 
cleation. The last is certainly the important characteristic but it has no direct 
connection with the specific phenomena of the multicomponent kinetics, which 
is the central point discussed here. This publication has to stimulate the exper- 
imental investigations in this direction because now it is absolutely clear what 
effects and characteristics are worth searching. 

In the last section the description of the "long nucleation" was done only for 
the period of nucleation and for the two component case. As for the descrip- 
tion of the second period it can be done in the same manner as for the case 
of the ordinary nucleation. It is necessary only to consider the two channels 
of nucleation in all places where we speak about the initial stage of condensa- 
tion and about the asymptotic stage of condensation. The situation when the 
evolution of the droplets from the first channel already attains the asymptotic 
stage (that's why we prefer to speak about stages of condensation instead of the 
periods of condensation) and the droplets 74 from the second channel are still 
at the initial stage is rather specific. But still in frames of approximations of 
initial stage and final stage it is easy to give the adequate description. 

For different components the power of consumption of the surplus vapor can 
be different. That's why at the given moment of time the different components 
can be at different stages. The evolution becomes unclear. In this context it is 
important to see the property of the global avalanche consumption, which states 
that the intermediate period has relatively (in comparison with the duration of 
the first period) short duration. Then for every concrete complex system one can 
act in frames of simple solutions and get the simple true kinetics of nucleation. 

Here we shall briefly summarize the new items proposed in this paper: 

• The consideration of kinetic and thermodynamic channels is given. 

• The correct description of the relaxation to the stationary state in the 
channel is given. 

• Essential errors of the elementary description of nucleation from [5] are 
corrected. 

• The description of the condensation of the main quantity of substance 
is given. This description corrects the description proposed in 9 . It is 
shown that the precise solution is absent. Several effective approximations 
are suggested. They are rather effective. 

• The situation of the "long nucleation" is analyzed. The methods to give 
the adequate description of this situation are proposed. 

But the main result of the article is the complete analytical description of the 
condensation process in the system with a mixture of substances in a metastable 

74 We have to consider the initial, intermediate and the final periods with reference for every 
component. 
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phase. All main characteristics of the process such as the numbers of droplets 
of dufferent sorts and the times of duration of nucleation in different channels 
of nucleation are determined explicitly. 

We have not given the estimates for the beginning of the transformation of 
the condensation process into the process of coalescense (super-condensation). 
The necessary estimate comes from condition that the coordinate of monodis- 
perse spectrum has to be 2 or 3 times greater than the critical size Vi c in the 
every corresponding channel. The announced estimate can be easily checked on 
the base of obtained solution, which gives the estimate for the time of the end 
of the condensation process 75 . 

The generalization of the description of the "long nucleation" for the mul- 
ticomponent case is rather simple. It is clear that in the n-component mixture 
condensation the maximum number of separate nucleation processes is limited 
by n. So, we have to act in the following manner: Consider the first channel 
of nucleation, then see what components aren't exhausted (the corresponding 
concentrations in kinetic channels are small). If there are some unexhausted 
components we can search for the channels where only these unexhausted com- 
ponents are the leading ones (the concentrations of all other components at 
thermodynamic channel are small). We shall do it in the sequence of increasing 
heights of activation barriers of nucleation. May be we shall find the channel 
with a necessary properties. This channel may not be the second in the total 
sequence of the activation barriers heights, but all channels before this channel 
are now out of action. The nucleation through this channel will be the sec- 
ondary nucleation and it is already described. Then we shall reconsider the set 
of unexhausted components (these components are the components, which are 
small in both kinetic channels) and seek for the kinetic channel with essential 
concentrations (concentrations, which aren't small) from the set of previously 
unexhausted components. This channel will be the leading one. This completes 
the loop, which can be repeated again. The number of such loops is limited 
from above by I. 

As the result we have the set of channels, which will be called as essential 
channels (both thermodynamic and kinetic). The nucleation in all channels 
from this set can be described analogously to the two component case. There 
is only one difficulty, which is specific to the multi component case. In the 
two component case the consumption of vapor by the droplets from the first 
channel occurs in the pseudo one-dimension manner (the concentration of one 
component is small, so the concentration of the other component is known also 
and, thus, all concentrations are known; moreover the value of the droplets 
volume is determined only by one component). In the multicomponent case 
there will be a real problem to describe this evolution. But with the help of 
methods presented in the section devoted to the second period this situation 
can be described analytically also. 

To give the necessary description one can split the evolution into the initial 

75 Also the statement from [5| that the theory of multicomponent coalescence is absent is 
wrong - the results presented in 1231 cover the case of multicomponent coalescence. 
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stage and the asymptotic stage. Since for droplets appeared from one channel 
the evolution is already at the final asymptotic stage and for the droplets from 
the other channel the evolution is at the initial stage at one and the same 
moment then it is preferable to speak about stages instead of periods. The 
initial stage and the final stage practically cover all evolution. This is the crucial 
point of description. Then it is rather easy to take into account the influence of 
the already existing spectrums of droplets appeared through the already closed 
channels, because we know the functional expressions for them as functions of 
time 76 . So, the functional form of equations will be the same as for the binary 
system and it is easy to get the approximate solution in analytical form. 

One have also to keep in mind that the final state corresponding to the differ- 
ent channels can be different. Certainly it doesn't produce essential difficulties. 
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